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This article reviews current efforts and plans for gravitational-wave detection, the gravitational- 
wave sources that might be detected, and the information that the detectors might extract from 
the observed waves. Special attention is paid to (i) the LIGO/ VIRGO network of earth-based, 
kilometer-scale laser interferometers, which is now under construction and will operate in the high- 
frequency band (1 to 10 4 Hz), and (ii) a proposed 5-million-kilometer-long Laser Interferometer 
Space Antenna (LISA), which would fly in heliocentric orbit and operate in the low- frequency 
band (1CT 4 to 1 Hz). LISA would extend the LIGO/VIRGO studies of stellar-mass (M ~ 2 
to 300Mq) black holes into the domain of the massive black holes (M ~ 1000 to 10 8 M Q ) that 
inhabit galactic nuclei and quasars. 



§1 Introduction 

Q According to general relativity theory, compact concen- 
trations of energy (e.g., neutron stars and black holes) 
should warp spacetime strongly, and whenever such an 
energy concentration changes shape, it should create a 
dynamically changing spacetime warpage that propagates 
out through the Universe at the speed of light. This prop- 
agating warpage is called a gravitational wave — a name 
that arises from general relativity's description of gravity 
as a consequence of spacetime warpage. 

Although gravitational waves have not yet been de- 
tected directly, their indirect influence has been seen and 
measured with such remarkable accuracy that their reality 
has been blessed even by the Nobel Prize Committee (that 
bastion of conservatism which explicitly denied Einstein 
the Prize for his relativity theories jlj): 

The 1993 Prize was awarded to Russell Hulse and 
Joseph Taylor for their discovery of the binary pulsar PSR 
1913+16 Q and for Taylor's observational demonstra- 
tion (with colleagues) S that the binary's two neutron 
stars are spiraling together at just the rate predicted by 
general relativity's theory of gravitational radiation reac- 
tion: from the observed orbit, one can compute the rate 
at which orbital energy should be emitted into gravita- 
tional radiation, and from this rate of energy loss one can 
compute the rate of orbital inspiral. The computed and 
observed inspiral rates agree to within the experimental 
accuracy, better than one per cent. 

Although this is a great triumph for Einstein, it is not a 
firm proof that general relativity is correct in all respects. 



1 This paper will be published in the Proceedings of the Snow- 
mass 95 Summer Study on Particle and Nuclear Astrophysics and 
Cosmology, eds. E. W. Kolb and R. Peccei (World Scientific, 
Singapore). 



Other relativistic theories of gravity (theories compatible 
with special relativity) predict the existence of gravita- 
tional waves; and some other theories predict the same 
inspiral rate for PSR 1913+16 as general relativity, to 
within the experimental accuracy ^, ||. Nevertheless, the 
experimental evidence for general relativity is so strong 
Q, that I shall assume it to be correct throughout this 
lecture except for occasional side remarks. 

There are a number of efforts, worldwide, to detect 
gravitational radiation. These efforts are driven in part 
by the desire to "see gravitational waves in the flesh," 
but more importantly by the goal of using the waves as a 
probe of the Universe and of the nature of gravity. And a 
powerful probe they should be, since they carry detailed 
information about gravity and their sources. 

There is an enormous difference between gravitational 
waves, and the electromagnetic waves on which our 
present knowledge of the Universe is based: 

• Electromagnetic waves are oscillations of the elec- 
tromagnetic field that propagate through spacetime; 
gravitational waves are oscillations of the "fabric" of 
spacetime itself. 

• Astronomical electromagnetic waves are almost al- 
ways incoherent superpositions of emission from in- 
dividual electrons, atoms, or molecules. Cosmic grav- 
itational waves are produced by coherent, bulk mo- 
tions of huge amounts of mass-energy — either mate- 
rial mass, or the energy of vibrating, nonlinear space- 
time curvature. 

• Since the wavelengths of electromagnetic waves are 
small compared to their sources (gas clouds, stellar 
atmospheres, accretion disks, ...), from the waves we 
can make pictures of the sources. The wavelengths 



of cosmic gravitational waves are comparable to or 
larger than their coherent, bulk-moving sources, so 
we cannot make pictures from them. Instead, the 
gravitational waves are like sound; they carry, in two 
independent waveforms, a stereophonic, symphony- 
like description of their sources. 

• Electromagnetic waves are easily absorbed, scattered, 
and dispersed by matter. Gravitational waves travel 
nearly unscathed through all forms and amounts of 
intervening matter |], [tJ . 

• Astronomical electromagnetic waves have frequencies 
that begin at / ~ 10 7 Hz and extend on upward by 
roughly 20 orders of magnitude. Astronomical grav- 
itational waves should begin at ~ 10 4 Hz (1000-fold 
lower than the lowest-frequency astronomical electro- 
magnetic waves), and should extend on downward 
from there by roughly 20 orders of magnitude. 

These enormous differences make it likely that: 

• The information brought to us by gravitational waves 
will be very different from (almost "orthogonal to" ) 
that carried by electromagnetic waves; gravitational 
waves will show us details of the bulk motion of 
dense concentrations of energy, whereas electromag- 
netic waves show us the thermodynamic state of op- 
tically thin concentrations of matter. 

• Most (but not all) gravitational-wave sources that 
our instruments detect will not be seen electro- 
magnetically, and conversely, most objects observed 
electromagnetically will never be seen gravitation- 
ally. Typical electromagnetic sources are stellar at- 
mospheres, accretion disks, and clouds of interstel- 
lar gas — none of which emit significant gravitational 
waves, while typical gravitational-wave sources are 
the cores of supernovae (which are hidden from elec- 
tromagnetic view by dense layers of surrounding stel- 
lar gas), and colliding black holes (which emit no elec- 
tromagnetic waves at all). 

• Gravitational waves may bring us great surprises. In 
the past, when a radically new window has been 
opened onto the Universe, the resulting surprises 
have had a profound, indeed revolutionary, impact. 
For example, the radio universe, as discovered in the 
1940s, 50s and 60s, turned out to be far more vio- 
lent than the optical universe; radio waves brought 
us quasars, pulsars, and the cosmic microwave radi- 
ation, and with them our first direct observational 
evidence for black holes, neutron stars, and the heat 
of the big bang ||, It is reasonable to hope that 
gravitational waves will bring a similar "revolution" . 



In this lecture I shall review the present status of at- 
tempts to detect gravitational radiation and plans for the 
future, and I shall describe some examples of information 
that we expect to garner from the observed waves. I shall 
begin, in Section ^|, with an overview of all the frequency 
bands in which astrophysical gravitational waves are ex- 
pected to be strong, the expected sources in each band, 
and the detection techniques being used in each. Then 
in subsequent sections I shall focus on (i) the "high fre- 
quency band" which is populated by waves from stellar 
mass black holes and neutron stars and is being probed 
by ground-based instruments: laser interferometers and 
resonant-mass antennas (Sections [| [| and ||), and (ii) 
the "low-frequency band" which is populated by waves 
from supermassive black holes and binary stars and is 
probed by space-based instruments: radio and optical 
tracking of spacecraft (Sections [?] and ||) . Finally, in Sec- 
tion H I shall describe the stochastic background of grav- 
itational waves that is thought to have been produced by 
various processes in the early universe, and prospects for 
detecting it in the various frequency bands. 

§2 Frequency Bands, Sources, and De- 
tection Methods 

Four gravitational-wave frequency bands are being ex- 
plored experimentally: the high-frequency band (HF; 
/ ~ 10 4 to 1 Hz), the low-frequency band (LF; / ~ 1 to 
10~ 4 Hz), the very-low frequency band (VLF; / ~ 10~ 7 
to 10~ 9 Hz), and the extremely-low-frequency band (ELF; 
/ - 10~ 15 to 10" 18 Hz). 

2.1 High- Frequency Band, 1 to 10 4 Hz 

A gravitational-wave source of mass M cannot be much 
smaller than its gravitational radius, 2GM/c 2 , and can- 
not emit strongly at periods much smaller than the light- 
travel time AttGM/c 3 around this gravitational radius. 
Correspondingly, the frequencies at which it emits are 
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where M is the mass of the Sun and G and c are New- 
ton's gravitation constant and the speed of light. To 
achieve a size of order its gravitational radius and thereby 
emit near this maximum frequency, an object presumably 
must be heavier than the Chandrasekhar limit, about 
the mass of the sun, Mq. Thus, the highest frequency 
expected for strong gravitational waves is / max ~ 10 4 
Hz. This defines the upper edge of the high-frequency 
gravitational-wave band. 

The high-frequency band is the domain of Earth-based 
gravitational-wave detectors: laser interferometers and 



resonant mass antennas. At frequencies below about 1 Hz, 
Earth-based detectors face nearly insurmountable noise 
(i) from fluctuating Newtonian gravity gradients (due, 
e.g., to the gravitational pulls of inhomogeneities in the 
Earth's atmosphere which move overhead with the wind), 
and (ii) from Earth vibrations (which are extremely diffi- 
cult to filter out mechanically below ~ f Hz). This defines 
the 1 Hz lower edge of the high-frequency band; to detect 
waves below this frequency, one must fly one's detectors 
in space. 

A number of interesting gravitational-wave sources fall 
in the high-frequency band: the stellar collapse to a neu- 
tron star or black hole in our Galaxy and distant galaxies, 
which sometimes triggers supernovae; the rotation and vi- 
bration of neutron stars (pulsars) in our Galaxy; the co- 
alescence of neutron-star and stellar-mass black-hole bi- 
naries (M & lOOOM©) in distant galaxies; and possibly 
such sources of stochastic background as vibrating loops 
of cosmic string, phase transitions in the early Universe, 
and the big bang in which the Universe was born. 

I shall discuss the high-frequency band in detail in Sec- 
tions 0- 



2.2 Low-Frequency Band, 10 4 to 1 Hz 

The low-frequency band, 10 -4 to f Hz, is the domain 
of detectors flown in space (in Earth orbit or in inter- 
planetary orbit). The most important of these are the 
Doppler tracking of spacecraft via microwave signals sent 
from Earth to the spacecraft and there transponded back 
to Earth (a technique that NASA has pursued since the 
early 1970's), and optical tracking of spacecraft by each 
other (laser interferometry in space, a technique now un- 
der development for possible flight in ~ 2014 or sooner). 

The 1 Hz upper edge of the low-frequency band is de- 
fined by the gravity-gradient and seismic cutoffs on Earth- 
based instruments; the ~ I0~ 4 Hz lower edge is defined 
by expected severe difficulties at lower frequencies in iso- 
lating spacecraft from the buffeting forces of fluctuating 
solar radiation pressure, solar wind, and cosmic rays. 

The low-frequency band should be populated by waves 
from short-period binary stars in our own Galaxy (main- 
sequence binaries, cataclysmic variables, white-dwarf bi- 
naries, neutron-star binaries, ...); from white dwarfs, neu- 
tron stars, and small black holes spiraling into massive 
black holes (M ~ 3 x f 5 to 3 x f 7 M Q ) in distant galax- 
ies; and from the inspiral and coalescence of supermassive 
black-hole binaries (M ~ fOO to I0 8 M o ). The upper 
limit, ~ IO 8 Af , on the masses of black holes that can 
emit in the low-frequency band is set by Eq. (|l|) with 
/ ib f0~ 4 Hz. There should also be a low-frequency 
stochastic background from such early-universe processes 
as vibrating cosmic strings, phase transitions, and the big- 
bang itself. 



I shall discuss the low-frequency band in detail in Sec- 
tions |?] and ||. 



2.3 Very-Low-Frequency Band, 10 7 to 10 9 Hz 

Joseph Taylor and others have achieved a remarkable 
gravity-wave sensitivity in the very-low-frequency band 
(VLF) by the timing of millisecond pulsars: When a gravi- 
tational wave passes over the Earth, it perturbs our rate of 
flow of time and thence the ticking rates of our clocks rel- 
ative to clocks outside the wave. Such perturbations will 
show up as apparent fluctuations in the times of arrival 
of the pulsar's pulses. If no fluctuations are seen at some 
level, we can be rather sure that neither Earth nor the 
pulsar is being bathed by gravitational waves of the cor- 
responding strength. If fluctuations with the same time 
evolution are seen simultaneously in the timing of several 
different pulsars, then the cause could well be gravita- 
tional waves bathing the Earth. 

By averaging the pulses' times of arrival over long peri- 
ods of time (months to tens of years) , a very high timing 
precision can be achieved, and correspondingly tight lim- 
its can be placed on the waves bathing the Earth or the 
pulsar. The upper edge of the VLF band, ~ I0~ 7 Hz, is 
set by the averaging time, a few months, needed to build 
up high accuracy; the lower edge, ~ 10 -9 Hz, is set by 
the time, ~ 20 years, since very steady millisecond pul- 
sars were first discovered. 



As we shall see in Section ^2 , strong gravitational- wave 
sources are generally compact, not much larger than their 
own gravitational radii. The only compact bodies that can 
radiate in the VLF band or below, i.e., at / I0~ 7 Hz, 
arc those with M £ 10 U M S [cf. Eq. (Q)]. Conventional 
astronomical wisdom suggests that compact bodies this 
massive do not exist, and that therefore the only strong 
waves in the VLF band and below are a stochastic back- 
ground produced by the same early-universe processes as 
might radiate at low and high frequencies: cosmic strings, 
phase transitions, and the big bang. 

Of course, conventional wisdom could be wrong. Nev- 
ertheless, it is conventional to quote measurement accu- 
racies in the VLF band and below in the language of a 
stochastic background: the fraction £l g (f) of the energy 
required to close the universe that lies in a bandwidth 
A/ = / centered on frequency /. The current 95%- 
confidence limit on Q g from pulsar timing in the VLF 
band is fi s (4 x 10~ 9 Hz) < 6 x 1Q- % H~ 2 where H is the 
Hubble constant in units of 100 km sec -1 Mpc" 1 Q. This 
is a sufficiently tight limit that it is beginning to cast 
doubt on the (not terribly popular) suggestion, that the 
Universe contains enough vibrating loops of cosmic string 
for their gravitational pulls to have seeded galaxy forma- 
tion coi, in. 
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Figure 1: The lines of force associated with the two polarizations of 
a gravitational wave. (From Ref. [114].) 
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Gravitational waves in the extremely-low-frequency band 
(ELF), 10~ 15 to 10 -18 Hz, should produce anisotropies in 
the cosmic microwave background radiation. The tightest 
limit from microwave observations comes from the lower 
edge of the ELF band / ~ 10 -18 Hz, where the gravita- 
tional wavelength is about tt times the Hubble distance, 
and the waves, by squeezing all of the space inside our cos- 
mological horizon in one direction, and stretching it all in 
another, should produce a quadrupolar anisotropy in the 
microwave background. The quadrupolar anisotropy mea- 
sured by the COBE satellite, if due primarily to gravita- 
tional waves (which it could be O), corresponds to an 
energy density fi g (Kr 18 Hz) ~ 10"T In Section 1 1 shall 
discuss the implications of this impressive ELF limit for 
the strength of the early-Universe stochastic background 
in the HF and LF bands. 

2.5 Other Frequency Bands and Other Detec- 
tion Methods 

A large number of other methods have been conceived of, 
for searching for gravitational radiation. Some of them 
would operate best in the HF, LF, VLF, and ELF bands 
described above; others would operate best at other fre- 
quencies. However, none has shown anywhere near the 
promise or the achievements of the methods described 
above (laser interferometry on Earth and in space, reso- 
nant mass antennas, Doppler tracking of spacecraft, tim- 
ing of pulsars, and anisotropy of microwave background) . 
For some references to other methods, see, e.g., [pi. 

§3 Ground-Based Laser Interferome- 
ters 

3.1 Wave Polarizations, Waveforms, and How 
an Interferometer Works 

According to general relativity theory (which I shall as- 
sume to be correct in this paper), a gravitational wave has 
two linear polarizations, conventionally called 4- (plus) 
and x (cross). Associated with each polarization there 
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Figure 2: Schematic diagram of a laser interferometer gravitational 
wave detector. (From Ref. ]14| . ) 



is a gravitational- wave field, /i+ or h x , which oscillates 
in time and propagates with the speed of light. Each 
wave field produces tidal forces (stretching and squeezing 
forces) on any object or detector through which it passes. 
If the object is small compared to the waves' wavelength 
(as is the case for ground-based interferometers and reso- 
nant mass antennas), then relative to the object's center, 
the forces have the quadrupolar patterns shown in Fig- 
ure [l]. The names "plus" and "cross" are derived from 
the orientations of the axes that characterize the force 
patterns [0]. 

A laser interferometer gravitational wave detector ( "in- 
terferometer" for short) consists of four masses that hang 
from vibration-isolated supports as shown in Figure |[ 
and the indicated optical system for monitoring the sepa- 
rations between the masses || Q . Two masses are near 
each other, at the corner of an "L", and one mass is at 
the end of each of the L's long arms. The arm lengths 
are nearly equal, L\ ~ L2 = L. When a gravitational 
wave, with frequencies high compared to the masses' ~ 1 
Hz pendulum frequency, passes through the detector, it 
pushes the masses back and forth relative to each other as 
though they were free from their suspension wires, thereby 
changing the arm-length difference, AL = L\ — L2. That 
change is monitored by laser interferometry in such a way 
that the variations in the output of the photodiode (the in- 
terferometer's output) are directly proportional to AL(t). 

If the waves are coming from overhead or underfoot and 
the axes of the + polarization coincide with the arms' 
directions, then it is the waves' + polarization that drives 
the masses, and AL(t)/L = h + (t). More generally, the 
interferometer's output is a linear combination of the two 
wave fields: 



AL{t) 
L 



= F+h+(t) + F x h x (t) = h(t) 



(2) 



The coefficients F + and F x are of order unity and depend 
in a quadrupolar manner on the direction to the source 
and the orientation of the detector |6) . The combination 
h(t) of the two /l's is called the gravitational-wave strain 



that acts on the detector; and the time evolutions of h(t), 
h+{t), and h x {t) are sometimes called waveforms. 

Interferometer test masses at present are made of trans- 
parent fused silica, though other materials might be used 
in the future. The masses' inner faces (shown white in 
Fig. ||) are covered with high-reflectivity dielectric coat- 
ings to form the indicated "mirrors", while the masses' 
outer faces are covered with anti-reflection coatings. The 
two mirrors facing each other on each arm form a Fabry- 
Perot cavity. A beam splitter splits a carefully prepared 
laser beam in two, and directs the resulting beams down 
the two arms. Each beam penetrates through the an- 
tireflection coating of its arm's corner mass, through the 
mass, and through the dielectric coating (the mirror); and 
thereby — with the length of the arm's Fabry-Perot cavity 
adjusted to be nearly an integral number of half wave- 
lengths of light — the beam gets trapped in the cavity. The 
cavity's end mirror has much higher reflectivity than its 
corner mirror, so the trapped light leaks back out through 
the corner mirror, and then hits the beam splitter where 
it recombines with light from the other arm. Most of 
the recombined light goes back toward the laser (where it 
can be returned to the interferometer by a "light-recycling 
mirror" labeled R), but a tiny portion goes toward the 
photodiode. 

When a gravitational wave hits the detector and moves 
the masses, thereby changing the lengths L\ and Li 
of the two cavities, it shifts each cavity's resonant fre- 
quency slightly relative to the laser frequency, and thereby 
changes the phase of the light in the cavity and the phase 
of the light that exits from the cavity toward the beam 
splitter. Correspondingly, the relative phase of the two 
beams returning to the splitter is altered by an amount 
A$ oc AL, and this relative phase shift causes a change 
in the intensity of the recombined light at the photodi- 
ode, A/pd oc A$ oc AL oc hit). Thus, the change of 
photodiode output current is directly proportional to the 
gravitational- wave strain h(t). This method of monitor- 
ing h(t), which was invented by Ronald Drever Jl5| as 
a modification of Rainer Weiss's Jl^] seminal concept for 
such an interferometer, is capable of very high sensitivity, 
as we shall see below. 

3.2 Wave Strengths and Interferometer Arm 
Lengths 

The strengths of the waves from a gravitational- 
wave source can be estimated using the "Newto- 
nian/quadrupole" approximation to the Einstein field 
equations. This approximation says that h ~ (G/c 4 )Q/r, 
where Q is the second time derivative of the source's 
quadrupole moment, r is the distance of the source from 
Earth (and G and c are Newton's gravitation constant and 
the speed of light) . The strongest sources will be highly 



nonspherical and thus will have Q ~ ML 2 , where M is 
their mass and L their size, and correspondingly will have 
Q ~ 2Mv 2 ~ where v is their internal velocity 

and E£? is the nonspherical part of their internal kinetic 
energy. This provides us with the estimate 



1 4G(££? n /c 2 ) 



(3) 



i.e., h is about 4 times the gravitational potential pro- 
duced at Earth by the mass-equivalent of the source's 
nonspherical, internal kinetic energy — made dimension- 
less by dividing by c 2 . Thus, in order to radiate strongly, 
the source must have a very large, nonspherical, internal 
kinetic energy. 

The best known way to achieve a huge internal ki- 
netic energy is via gravity; and by energy conservation 
(or the virial theorem), any gravitationally-induced ki- 
netic energy must be of order the source's gravitational 
potential energy. A huge potential energy, in turn, re- 
quires that the source be very compact, not much larger 
than its own gravitational radius. Thus, the strongest 
gravity-wave sources must be highly compact, dynamical 
concentrations of large amounts of mass (e.g., colliding 
and coalescing black holes and neutron stars). 

Such sources cannot remain highly dynamical for long; 
their motions will be stopped by energy loss to gravita- 
tional waves and/or the formation of an all-encompassing 
black hole. Thus, the strongest sources should be tran- 
sient. Moreover, they should be very rare — so rare that 
to see a reasonable event rate will require reaching out 
through a substantial fraction of the Universe. Thus, just 
as the strongest radio waves arriving at Earth tend to be 
extragalactic, so also the strongest gravitational waves are 
likely to be extragalactic. 

For highly compact, dynamical objects that radiate in 
the high-frequency band, e.g. colliding and coalescing neu- 
tron stars and stellar-mass black holes, the internal, non- 
spherical kinetic energy E£?/c 2 is of order the mass of 
the Sun; and, correspondingly, Eq. (||) gives h ~ 10~ 22 
for such sources at the Hubble distance (3000 Mpc, i.e., 
10 10 light years); h ~ 10" 21 at 200 Mpc (a best-guess dis- 
tance f or se veral neutron-star coalescences per year; see 
Section 5.2), h ~ 10~ 20 at the Virgo cluster of galaxies 
(15 Mpc); and h ~ 10~ 17 in the outer reaches of our own 
Milky Way galaxy (20 kpc). These numbers set the scale 
of sensitivities that ground-based interferometers seek to 
achieve: h ~ 10~ 21 to 10~ 22 . 

When one examines the technology of laser interfer- 
ometry, one sees good prospects to achieve measurement 
accuracies AL ~ 10 -16 cm (1/1000 the diameter of the 
nucleus of an atom) . With such an accuracy, an interfer- 
ometer must have an arm length L = AL/h ~ 1 to 10 km, 
in order to achieve the desired wave sensitivities, 10~ 21 to 



Figure 3: Artist's conception of one of the LIGO interferometers. 
[Courtesy the LIGO Project.] 

10~ 22 . This sets the scale of the interferometers that are 
now under construction. 

3.3 LIGO, VIRGO, and the International In- 
ter jerometric Network 

Interferometers are plagued by non-Gaussian noise, e.g. 
due to sudden strain releases in the wires that suspend 
the masses. This noise prevents a single interferometer, 
by itself, from detecting with confidence short-duration 
gravitational-wave bursts (though it might be possible for 
a single interferometer to search for the periodic waves 
from known pulsars). The non-Gaussian noise can be 
removed by cross correlating two, or preferably three 
or more, interferometers that are networked together at 
widely separated sites. 

The technology and techniques for such interferometers 
have been under development for nearly 25 years, and 
plans for km-scale interferometers have been developed 
over the past 14 years. An international network con- 
sisting of three km-scale interferometers, at three widely 
separated sites, is now in the early stages of construc- 
tion. It includes two sites of the American LIGO Project 
( "Laser Interferometer Gravitational Wave Observatory" ) 
PI , and one site of the French/Italian VIRGO Project 
(named after the Virgo cluster of galaxies) [Il7} . 

LIGO will consist of two vacuum facilities with 4- 
kilometer-long arms, one in Hanford, Washington (in the 
northwestern United States; Fig. ||) and the other in Liv- 
ingston, Louisiana (in the southeastern United States). 
These facilities are designed to house many successive gen- 
erations of interferometers without the necessity of any 
major facilities upgrade; and after a planned future ex- 
pansion, they will be able to house several interferometers 
at once, each with a different optical configuration opti- 
mized for a different type of wave (e.g., broad-band burst, 
or narrow-band periodic wave, or stochastic wave). The 
LIGO facilities and their first interferometers are being 



constructed by a team of about 80 physicists and engineers 
at Caltech and MIT, led by Barry Barish (the PI) and 
Gary Sanders (the Project Manager). Robbie Vogt (who 
directed the project during the pre-construction phase) is 
in charge of the final design and construction of LIGO's 
first interferometers, Stan Whitcomb is in charge of inter- 
ferometer R&D, and Albert Lazzarini is the system engi- 
neer and Rai Weiss the cognizant scientist for integration 
of all parts of LIGO. 

A number of other research groups are making impor- 
tant contributions to LIGO: Bob Byers' group at Stan- 
ford is developing Nd:YAG lasers, Peter Saulson's group 
at Syracuse and Vladimir Braginsky's group in Moscow 
are developing test-mass suspension systems and studying 
noise in them; Jim Faller's group at JILA is developing 
active vibration isolation systems; Ron Drever's group at 
Caltech is developing advanced interferometers; and Sam 
Finn's group at Northwestern and my group at Caltech 
are developing data analysis techniques. A number of 
other groups are likely to join the LIGO effort in the next 
few years. A formal association of LIGO-related scientists 
(the LIGO Research Community, an analog of a "user's 
group" in high-energy physics) is being organized, and a 
LIGO Program Advisory Committee will be formed soon, 
with voting membership restricted to people outside the 
Caltech/MIT LIGO team, to advise the LIGO manage- 
ment. 

The VIRGO Project is building one vacuum facility in 
Pisa, Italy, with 3-kilometer-long arms. This facility and 
its first interferometers are a collaboration of more than a 
hundred physicists and engineers at the INFN (Frascati, 
Napoli, Perugia, Pisa), LAL (Orsay), LAPP (Annecy), 
LOA (Palaiseau), IPN (Lyon), ESPCI (Paris), and the 
University of Illinois (Urbana), under the leadership of 
Alain Brillet and Adalberto Giazotto. 

Both LIGO and VIRGO are scheduled for completion in 
the late 1990s, and their first gravitational-wave searches 
are likely to be performed in 2000 or 2001. 

LIGO alone, with its two sites which have parallel arms, 
will be able to detect an incoming gravitational wave, 
measure one of its two waveforms, and (from the time 
delay between the two sites) locate its source to within 
a ~ 1° wide annulus on the sky. LIGO and VIRGO to- 
gether, operating as a coordinated international network, 
will be able to locate the source (via time delays plus the 
interferometers' beam patterns) to within a 2-dimensional 
error box with size between several tens of arcminutes and 
several degrees, depending on the source direction and on 
the amount of high-frequency structure in the waveforms. 
They will also be able to monitor both waveforms h + (t) 
and h x (t) (except for frequency components above about 
1kHz and below about 10 Hz, where the interferometers' 
noise becomes severe). 

The accuracies of the direction measurements and the 



ability to monitor more than one waveform will be severely 
compromised when the source lies anywhere near the 
plane formed by the three LIGO/ VIRGO interferome- 
ter locations. To get good all-sky coverage will require a 
fourth interferometer at a site far out of that plane; Japan 
and Australia would be excellent locations, and research 
groups there are carrying out research and development 
on interferometric detectors, aimed at such a possibility. 
A 300-meter prototype interferometer called TAMA is un- 
der construction in Tokyo, and a 400-meter prototype 
called AIGO400 has been proposed for construction north 
of Perth. 

Two other groups are major players in this field, one 
in Britain led by James Hough, the other in Germany, 
led by Karsten Danzmann. These groups each have 
two decades of experience with prototype interferometers 
(comparable experience to the LIGO team and far more 
than anyone else) and great expertise. Frustrated by in- 
adequate financing for a kilometer-scale interferometer, 
they are constructing, instead, a 600 meter system called 
GEO600 near Hannover, Germany. Their goal is to de- 
velop, from the outset, an interferometer with the sort 
of advanced design that LIGO and VIRGO will attempt 
only as a "second-generation" instrument, and thereby 
achieve sufficient sensitivity to be full partners in the 
international network's first gravitational-wave searches; 
they then would offer a variant of their interferometer as 
a candidate for second-generation operation in the much 
longer arms of LIGO and/or VIRGO. It is a seemingly 
audacious plan, but with their extensive experience and 
expertise, the British/German collaboration might pull it 
off successfully. 

3.4 Interferometer Development and Noise 
Sources 

It is not possible, in the short vacuum systems now avail- 
able (arm lengths < 40 meters), to develop and test a 
multikilomcter interferometer as a single unit. This is 
because the various noise sources that plague an inter- 
ferometer scale differently from each other with length 
L and with gravity- wave frequency /. As a result, the 
various components of the multikilomcter interferometers, 
and the various techniques to be used in them, are being 
developed and tested separately in a number of differ- 
ent laboratories, and will only be combined together into 
a single interferometer when the LIGO/VIRGO vacuum 
systems are completed. 

The best known of the LIGO-Project laboratories in 
which components and techniques are being developed is 
the 40-meter prototype interferometer at Caltech (Fig. 
^). This prototype focuses on the development of meth- 
ods and components to control "displacement noise," i.e., 
those noise sources that push the mirrored test masses 



Figure 4: The LIGO Project's 40-meter "Mark II" prototype inter- 
ferometer at Caltech. This prototype went into operation in 1993. 
It has much larger vacuum chambers, to accommodate bigger and 
better seismic isolation stacks, than those of the previous "Mark I" 
prototype (which operated from the early 1980s to 1992). [Courtesy 
the LIGO Project.] 

back and forth as would a gravity wave. The princi- 
pal sources of displacement noise are seismic vibrations 
of the ground beneath the interferometer (which are fil- 
tered out by the masses' suspension wires and by "isola- 
tion stacks" made of successive layers of steel and rub- 
ber), and thermally-induced vibrations of the test masses 
and of the wires that suspend them (vibrations that are 
controlled by designing the test masses and suspensions 
with great care and constructing them from low-loss, i.e. 
high "Q", materials). 

Among the LIGO Project's other laboratories, there is 
a shorter-armed prototype-interferometer facility at MIT, 
which is devoted to developing methods and components 
for controlling noise in the phase of the interferometer's 
light beams. Since the gravity wave makes itself known 
by the phase shift that it puts on the light of one in- 
terferometer arm relative to the other, this phase noise 
can simulate a gravity wave. Among the various causes 
of phase noise, the one that is the most fundamental is 
photon shot noise due to the random times at which the 
light's photons arrive at the photodiode (cf. Fig. ||). 

Once the myriad of other noise sources have been 
brought under control, shot noise, thermal noise (i.e., 
thermally induced vibrations), and seismic noise (i.e., 
ground vibrations) are likely to be the ultimate impedi- 
ments to detecting and studying gravitational waves. Fig- 
ure [5] shows the spectra expected for each of these three 
noises in the first interferometers that will operate in 
LIGO. At frequencies above 200 Hz, shot noise dominates; 
between 200 Hz and 40 Hz, thermal noise in the suspen- 
sion wires dominates; and below 40 Hz, seismic noise dom- 
inates. 

During LIGO's operations, step-by-step improvements 
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Figure 5: The expected noise spectrum in each of LIGO's first 4-km 
interferometers (upper solid curve) and in more advanced interfer- 
ometers (lower solid curve). The dashed curves show various con- 
tributions to the first interferometers' noise. Plotted horizontally is 
gravity wave frequency /; plotted vertically is h(f), the square root 
of the spectral density of the detector's output h(t) = AL(t)/L in 
the absence of a gravity wave. The rms h noise in a bandwidth A/ 
at frequency / is ft rms = h(f)y/ A/. (From Ref. |E3|.) 

will be made in the control of these three noise sources 
[0, thereby pushing the overall noise spectrum down- 
ward from the "first-interferometer" level toward the 
"advanced-interferometer" level shown in Figure |. As 
we shall see below, the sensitivity of the first interferom- 
eters might be inadequate to detect gravitational waves. 
However, we are quite confident that at some point during 
the improvement from first interferometers to advanced, 
a plethora of gravitational waves will be found and will 
start bringing us exciting information about fundamental 
physics and the Universe. 

Notice from Figure [| that the advanced LIGO inter- 
ferometers are expected to have their optimal sensitivity 
at / ~ 100 Hz, and rather good sensitivity all the way 
from / ^ 10 Hz at the low-frequency end to / ~ 500 Hz 
at the high-frequency end. Below 10 Hz, seismic noise, 
creeping through the isolation stacks, will overwhelm all 
gravitational-wave signals; and above 500 Hz, photon shot 
noise may overwhelm the signals. 

Figure || gives an impression of the present state of in- 
terferometer technology and the rate of progress. This 
figure shows a sequence of noise spectra in the 40-meter 
prototype interferometer at Caltech, during 1990-94 when 
the 40-meter R&D emphasis was on improving the low- 
frequency noise performance. The top two (noisiest) spec- 
tra are snapshots of the original "Mark I" prototype per- 
formance in October 1990 and June 1992; the lower two 
(quieter) spectra are from the rebuilt, "Mark II" pro- 
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Figure 6: Measured noise spectra in the Caltech 40-meter prototype 
interferometer (Fig. U). Since this prototype is devoted tq_ learning 
to control displacement noise, the spectra shown are AL(/), the 
square root of the spectral density of the measured arm-length dif- 
ference. Each of the many spectral lines is well understood, and 
most could be removed if their removal were of high priority (e.g., 
they are multiples of the 60 Hz line frequency sneaking into the 
apparatus by known routes). Those few, very narrow lines that 
cannot be physically removed by cleaning up the instrument (e.g., 
thermally-driven violin-mode resonances of the wires that suspend 
the test masses) will be filtered out in the data analysis. Thus, the 
interferometer sensitivity is characterized by theLXpntinuum noise 
floor and not the lines. (From Ref. faq].) 

totype of Figure [|, in March and October 1994. The 
smooth, solid line, for comparison, is the displacement 
noise goal for the first 4-km interferometers in LIGO (i.e., 
it is the upper solid curve of Fig. |^). 

Note that the prototypes's arm-length difference was 
being monitored, in October 1994, to within an rms 
noise level (in a bandwidth equal to frequency) AL rms = 
V / /AL(/) < 8 x 10~ 16 cm over the frequency range 200 to 
1000 Hz. This corresponds to an rms gravitational-wave 
noise level h rms = AL rms /40m < 2 x 10 -19 , the best that 
any gravitational-wave detector has yet achieved. 

3. 5 Semiquantitative Discussion of Interferom- 
eter Noise 

The LIGO and VIRGO interferometers are expected to 
have rms noise levels h tms 10~ 22 corresponding to test- 
mass position noises AL rms <i hL ~ 10~ 16 cm. 10~ 16 cm 
is awfully small: 1/1000 the diameter of the nucleus of an 
atom, and 10~ 12 the wavelength of the light being used to 
monitor the masses' motions. How can one possibly mon- 
itor such small motions? The following estimate explains 
how. 

One adjusts the reflectivities of the interferometer's cor- 
ner mirrors so the two arms store the laser light on average 
for about half a cycle of a ~ 100 Hz gravitational wave, 
which means for ~ 100 round trips. The light in each arm 



thereby acquires a phase shift 

A$ - 100 x 4ttAL/A - 10 -9 , (4) 

where A ~ 10 cm is the wavelength of light. If the inter- 
ference of the light from the two beams is done optimally, 
then this phase shift (equal and opposite in the two arms) 
can be measured at the photodiode to an accuracy that is 
governed by the light's photon shot noise, A<E> ~ l/vW, 
where N is the number of photons that enter the inter- 
ferometer from the laser during the half-cycle of photon 
storage time. (This 1/yN is the usual photon fluctuation 
in a quantum mechanical "coherent state" of light.) Thus, 
to achieve the required accuracy, A<& ~ 10~ 9 , in the face 
of photon shot noise, requires N ~ 10 18 photons in 0.01 
second, which means a laser power of ~ 100 Watts. 

By cleverness |l^], one can reduce the required laser 
power: The light is stored in the interferometer arms for 
only a half gravity-wave period (~ 100 round trips) be- 
cause during the next half period the waves would reverse 
the sign of AL, thereby reversing the sign of the phase 
shift being put onto the light and removing from the light 
the signal that had accumulated in the first half period. 
In just 100 round trips, however, the light is attenuated 
hardly at all. One therefore reuses the light, over and 
over again. This is done by (i) operating the interferome- 
ter with only a tiny fraction of the recombined light going 
out toward the photodiode, and almost all of it instead 
going back toward the laser, and by (ii) placing a mirror 
(marked R in Fig. ||) between the laser and the interfer- 
ometer in just such a position that the entire interferom- 
eter becomes an optical cavity driven by the laser — with 
its arms as two subcavities. Then the mirror R recycles 
the recombined light back into the interferometer in phase 
with the new laser light, thereby enabling a laser of, say, 
5 Watts to behave like one of 100 Watts or more. 

Turn from photon shot noise to thermal noise. How, 
one might ask, can one possibly expect to monitor 
the mirrors' motions at a level of 10 -16 cm when the 
room-temperature atoms of which the fused-silica mir- 
rors are made vibrate thermally with amplitudes Ai rms = 
^/kT/muj 2 ~ 10~ 10 cm? (Here k is Boltzmann's constant, 
T is room temperature, m is the atomic rest mass, and 
u ~ 10 14 s _1 is the angular frequency of atomic vibra- 
tion.) The answer is that these individual atomic vibra- 
tions are unimportant. The light beam, with its ~ 5cm 
diameter, averages over the positions of ~ 10 17 atoms 
in the mirror, and with its 0.01s storage time it aver- 
ages over ~ 10 11 vibrations of each atom. This spatial 
and temporal averaging makes the vibrations of individ- 
ual atoms irrelevant. Not so irrelevant, however, are the 
lowest-frequency normal-mode vibrations of the mirror- 
endowed masses (since these modes experience much less 
time averaging than the faster atomic vibrations). Assum- 
ing a mass m ~ (a few tens of kg), these normal modes 



Figure 7: The seismic isolation stack that was recently installed in 
the LIGO Project's Mark II prototype interferometer at Caltech. 
When the interferometer is in operation, a small tower is mounted 
on the top steel plate and from the tower hangs one of the interfer- 
ometer's mirror-endowed masses. [Courtesy the LIGO Project.] 

have angular frequencies u> ~ 10 s , so their rms vibra- 
tion amplitude is A/ rms = ^/kT/muj 2 ~ 10~ 14 cm. This 
is 100 times larger than the signals we wish to monitor; 
but if these modes have high quality factors (high Q's; low 
losses), then the vibrations will be very steady over the 
interferometer's averaging time of 0.01 s, and correspond- 
ingly, their effects will average down by more than a factor 
100. Similar considerations apply to the thermal noise in 
the masses' suspension wires. For detailed discussions of 
fascinating and not-fully-understood physics issues that 
influence the thermal noise, see, e.g., Refs. 0, ^o|, 53. 

At the LIGO sites, and most any other quiet location 
on Earth, the ground is continually shaking with an rms 
displacement Al lms ~ 10~ 8 cm (100 Hz//) 3 / 2 . This is 
10 7 times larger than the motions one seeks to monitor. 
At frequencies above 10 Hz or so, one can protect the 
masses from these seismic vibrations by simple (but care- 
fully designed) passive isolation stacks. Each element in 
the stack is a mass and a spring (a harmonic oscillator) 
with normal- mode frequency fo ~ (a few Hz). When seis- 
mic noise tries to drive this harmonic oscillator far above 
its resonant frequency [in our case at / <; (a few tens of 
Hz)], the amplitude of its response is attenuated relative 
to the driving motion by a factor (/o//) 2 [in our case a 
factor <; 10 2 ]. Thus, each oscillator in the stack will pro- 
vide a reduction <; 10 2 in AZ rms , so a stack of four or five 
oscillators is enough to provide the required isolation. Fig- 
ure ^ shows an isolation stack — made of four steel plates 
and four sets of viton rubber springs (not quite visible 
between the plates) — that is now operating in the Mark 
II prototype interferometer of Figure [|. This stack and 
the pendulum wires that suspend the mirror-endowed test 
masses provide five layers of isolation. The installation of 
this new stack was responsible for the sharp drop in low- 
frequency noise in Figure |^ between June 1992 (Mark I) 



and March 1994 (Mark II). 

The above rough estimates suggest (as Weiss realized as 
early as 1972 jl(|) that it is possible for interferometers 
to achieve the required sensitivities, h lms ~ 10~ 22 and 
AL ~ 10 -16 cm. However, going from these rough esti- 
mates to a real working interferometer, and doing so in the 
face of a plethora of other noise sources, is a tremendous 
experimental challenge — one that has occupied a number 
of excellent experimental physicists since 1972. 

§4 Resonant-Mass Antennas 

A resonant-mass antenna for gravitational radiation con- 
sists of a solid body that (heuristically speaking) rings 
like a bell when a gravitational wave hits it. This body 
(the resonant mass) is usually a cylinder, but future vari- 
ants are likely to be spheres or sphere-like, e.g. a trun- 
cated icosahedron gravitational-wave antenna or TIGA 
PH - The resonant mass is typically made from an al- 
loy of aluminum and weighs several tons, but some have 
been made of niobium or single-crystal silicon or sapphire 
(but with masses well below a ton). To control thermal 
noise, the resonant mass is usually cooled cryogenically to 
liquid-helium temperatures or below. 

The resonant-mass antenna is instrumented with an 
electromagnetic transducer and electronics, which mon- 
itor the complex amplitude of one or more of the mass's 
normal modes. When a gravitational wave passes through 
the mass, its frequency components near each normal- 
mode frequency f a drive that mode, changing its com- 
plex amplitude; and the time evolution of the changes 
is measured within some bandwidth A/ by the trans- 
ducer and electronics. Current resonant-mass antennas 
are narrow-band devices (A/// Q <C 1) but in the era of 
LIGO/VIRGO, they might achieve bandwidths as large 
as Af/f ~ 1/3. 

Resonant-mass antennas for gravitational radiation 
were pioneered by Joseph Weber about 35 years ago p| , 
and have been pushed to ever higher sensitivity by Weber 
and a number of other research groups since then. For ref- 
erences and an overview of the present and future of such 
antennas see, e.g., Ref. At present there is a network 
of such antennas [^5| , cooled to 3K, and operating with 
an rms noise level for broad-band gravity-wave bursts of 
h rms ~ 6 x 10~ 19 . The network includes an aluminum 
cylinder called EXPLORER built by a group at the Uni- 
versity of Rome, Italy (Edoardo Amaldi, Guido Pizclla, 
et. al.); an aluminum cylinder at Louisiana State Univer- 
sity, USA (Bill Hamilton, Warren Johnson, et. al.); and 
a niobium cylinder at the University of Perth, Australia 
(David Blair et. al.). This network has been in operation, 
searching for waves, for several years. 

The next generation of resonant-mass antennas is now 
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Figure 8: The rms noise curves h(f) (measured in strain per root 
Hz) for a "xylaphone" of TIGA gravitational- wave detectors (solid 
curves) for signals of random polarization and direction p2| . The 
TIGA's are presumed instrumented and cooled sufficiently well that 
their sensitivity is at the standard quantum limit. Their central 
frequencies, radii and masses (assuming aluminum material) are {1.0 
kHz, 1.30 m, 25.1 ton}, {1.25 kHz, 1.04 m, 12.8 ton}, {1.50 kHz, 
0.87 m, 7.4 ton}, {1.75 kHz, 0.74 m, 4.7 ton}, {2.0 kHz, 0.65 m, 
3.1 ton}. Shown for comparison are the noise curves for the first 
LIGO interferometer with random wave polarization and direction 
multiplied by y/E (dotted curve; Fig. ^|), and for the first GFO600 
detector operated in a narrow-band mode (dashed curve; Ref. Pa]). 

under construction at the University of Rome (NAU- 
TILUS) and at the University of Legarno, Italy (AU- 
RIGA). These are several-ton aluminum bars cooled to 
0.05K; their rms design sensitivities for wave bursts are 
(several) x 10" 20 ||. 

A subsequent generation, which hopefully would op- 
erate in the LIGO/VIRGO era, is being discussed and 
planned Q. These are 1 to 100 ton spheres or TIGA's 
cooled to ~ 0.01 — 0.05K, with sensitivity goals of ~ 
10~ 21 . Such antennas might be built by an American col- 
laboration, a Brazilian collaboration, an Italian collabora- 
tion called "Omega", and/or a Dutch collaboration called 
"Grail" . Their spherical or near-spherical shapes make 
them omnidirectional and should give them several-times 
higher sensitivities than can be achieved by cylinders at 
the same frequency. 

The attractiveness of such antennas in the LIGO/ 
VIRGO era lies in their ability to operate with impressive 
sensitivity in the uppermost reaches of the high-frequency 
band, ~ 10 3 to 10 4 Hz, where photon shot noise debilitates 
the performance of interferometric detectors (cf. Fig. ||). 
Figure || shows the projected rms noise curves of a family 
of TIGA detectors, each instrumented to operate at the 
"standard quantum limit" for such a detector (a nontriv- 
ial experimental task). Shown for comparison is the rms 
noise of the first LIGO interferometer — which, of course, 
is not optimized for the kHz band. The GEO600 inter- 
ferometer, with its advanced design, can be operated in 
a narrow-band, high-frequency mode (and probably will 
be so operated in ~ 1999. Its rms design sensitivity in 
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such a mode is also shown in Figure g. The TIGA sen- 
sitivities are sufficiently good in the kHz band, compared 
to early L1GO and GEO interferometers, that, although 
they probably cannot begin to operate until somewhat af- 
ter the beginning of the LIGO /VIRGO era, they might be 
fully competitive when they do operate, and might play 
an important role in the kHz band. 

§5 High-Frequency Gravitational- Wave 
Sources: Coalescing Compact Binaries 

The best understood of all gravitational-wave sources are 
coalescing, compact binaries composed of neutron stars 
(NS) and black holes (BH). These NS/NS, NS/BH, and 
BH/BH binaries may well become the "bread and butter" 
of the LIGO/ VIRGO diet. 

The Hulse-Taylor f| § binary pulsar, PSR 1913+16, 
is an example of a NS/NS binary whose waves could 
be measured by LIGO/VIRGO, if we were to wait long 
enough. At present PSR 1913+16 has an orbital fre- 
quency of about 1/(8 hours) and emits its waves predom- 
inantly at twice this frequency, roughly 1CV 4 Hz, which is 
in the low-frequency band — far too low to be detected by 
LIGO/VIRGO. However, as a result of their loss of orbital 
energy to gravitational waves, the PSR 1913+16 NS's are 
gradually spiraling inward. If we wait roughly 10 8 years, 
this inspiral will bring the waves into the LIGO/VIRGO 
high-frequency band. As the NS's continue their inspiral, 
the waves will then sweep upward in frequency, over a time 
of about 15 minutes, from 10 Hz to ~ 10 3 Hz, at which 
point the NS's will collide and coalesce. It is this last 
15 minutes of inspiral, with ~ 16, 000 cycles of waveform 
oscillation, and the final coalescence, that LIGO/VIRGO 
seeks to monitor. 

5.1 Wave Strengths Compared to LIGO Sensi- 
tivities 

Figure ^ compares the projected sensitivities of inter- 
ferometers in LIGO 0| with the wave strengths from 
the last few minutes of inspiral of BH/BH, NS/BH, and 
NS/NS binaries at various distances from Earth. The two 
solid curves at the bottoms of the stippled regions (labeled 
h-rms) are the rms noise levels for broad-band waves that 
have optimal direction and polarization. The tops of the 
stippled regions (labeled hsB for "sensitivity to bursts") 
are the sensitivities for highly confident detection of ran- 
domly polarized, broad-band waves from random direc- 
tions (i.e., the sensitivities for high confidence that any 
such observed signal is not a false alarm due to Gaus- 
sian noise). The upper stippled region and its bounding 
curves are the expected performances of the first interfer- 
ometers in LIGO; the lower stippled region and curves are 
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Figure 9: LIGO's projected broad-band noise /i rms and sensitivity 
to bursts hsB (Fig. ^ and Ref. Jl4| ) compared with the strengths of 
the waves from the last few minutes of inspiral of compact binaries. 
The signal to noise ratios are y2 higher than in Ref. ]14| because of 
a factor 2 error in Eq. (29) of Ref. 



performances of more advanced LIGO interferometers; cf. 
Figure ||. 

As the NS's and/or BH's spiral inward, their waves 
sweep upward in frequency (left to right in the dia- 
gram) . The dashed lines show their "characteristic" signal 
strength h c (approximately the amplitude h of the waves' 
oscillations multiplied by the square root of the number 
of cycles spent near a given frequency, \fn); the signal- 
to-noise ratio is this h c divided by the detector's V5h lms , 
S/N — h c /(\/~5h rrns ), where the y/b converts h ims from 
"optimal direction and polarization" to "random direction 
and polarization") |L4| ||. The arrows along the bottom 
inspiral track indicate the time until final coalescence for 
an NS/NS binary and the separation between the NS cen- 
ters of mass. Each NS is assumed to have a mass of 1.4 
suns and a radius ~ 10 km; each BH, 10 suns and ~ 20 
km. 

Notice that the signal strengths in Figure ^ are in good 
accord with our rough estimates based on Eq. (|^); at the 
endpoint (right end) of each inspiral, the number of cy- 
cles n spent near that frequency is of order unity, so the 
quantity plotted, h c ~ hy/n, is about equal to h — and at 
distance 200 Mpc is roughly 10 -21 , as we estimated in 
Section |J. 

5.2 Coalescence Rates 

Such final coalescences are few and far between in our 
own galaxy: about one every 100,000 years, according to 
1991 estimates by Phinney p^] and by Narayan, Piran, 
and Shemi p8[, based on the statistics of binary pulsar 
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searches in our galaxy which found three that will coa- 
lesce in less than 10 10 years. Extrapolating out through 
the universe on the basis of the density of production of 
blue light (the color produced predominantly by massive 
stars), Phinney |^7| and Narayan et. al. |^8| infer that to 
see several NS/NS coalescences per year, LIGO/VIRGO 
will have to look out to a distance of about 200 Mpc (give 
or take a factor ~ 2); cf. the "NS/NS inspiral, 200 Mpc" 
line in Figure |^. Since these estimates were made, the 
binary pulsar searches have been extended through a sig- 
nificantly larger volume of the galaxy than before, and no 
new ones with coalescence times <, 10 10 years have been 
found; as a result, the binary-pulsar-search-based best es- 
timate of the coalescence rate should be revised downward 
p9| , perhaps to as little as one every million years in our 
galaxy, corresponding to a distance 400 Mpc for several 
per year p9| . 

A rate of one every million years in our galaxy is ~ 1000 
times smaller than the birth rate of the NS/NS binaries' 
progenitors: massive, compact, main-sequence binaries 
p7| , |2g|| . Therefore, either 99.9 per cent of progenitors 
fail to make it to the NS/NS state (e.g., because of bi- 
nary disruption during a supernova or forming TZO's), 
or else they do make it, but they wind up as a class of 
NS/NS binaries that has not yet been discovered in any 
of the pulsar searches. Several experts on binary evolu- 
tion have argued for the latter @ fl], |||: most NS/NS 
binaries, they suggest, may form with such short orbital 
periods that their lifetimes to coalescence are significantly 
shorter than normal pulsar lifetimes (~ 10 7 years); and 
with such short lifetimes, they have been missed in pul- 
sar searches. By modeling the evolution of the galaxy's 
binary star population, the binary experts arrive at best 
estimates as high as 3 x 10 -4 coalescences per year in our 
galaxy, corresponding to several per year out to 60 Mpc 
distance |30) . Phinney J27| describes other plausible pop- 
ulations of NS /NS binaries that could increase the event 
rate, and he argues for "ultraconservative" lower and up- 
per limits of 23 Mpc and lOOOMpc for how far one must 
look to see several coalescence per year. 

By comparing these rate estimates with the signal 
strengths in Figure [| we see that (i) the first interfer- 
ometers in LIGO/VIRGO have a possibility but not high 
probability of seeing NS/NS coalescences; (ii) advanced 
interferometers are almost certain of seeing them (the 
requirement that this be so was one factor that forced 
the LIGO/VIRGO arm lengths to be so long, several 
kilometers); and (iii) they are most likely to be discov- 
ered roughly half-way between the first and advanced 
interferometers — which means by an improved variant of 
the first interferometers several years after LIGO opera- 
tions begin. 

We have no good observational handle on the coales- 
cence rate of NS/BH or BH/BH binaries. However, the- 



ory suggests that their progenitors might not disrupt dur- 
ing the stellar collapses that produce the NS's and BH's, 
so their coalescence rate could be about the same as the 
birth rate for their progenitors: ~ 1/100, 000 years in our 
galaxy. This suggests that within 200 Mpc distance there 
might be several NS/BH or BH/BH coalescences per year. 
p7[ p8[ [5(j| ^2). This estimate should be regarded as a 
plausible upper limit on the event rate and lower limit on 
the distance to look |27], |2^| . 

If this estimate is correct, then NS/BH and BH/BH 
binaries will be seen before NS/NS, and might be seen by 
the first LIGO/VIRGO interferometers or soon thereafter; 
cf. Figure ^. However, this estimate is far less certain than 
the (rather uncertain) NS/NS estimates! 

Once coalescence waves have been discovered, each fur- 
ther improvement of sensitivity by a factor 2 will increase 
the event rate by 2 3 ~ 10. Assuming a rate of several 
NS/NS per year at 200 Mpc, the advanced interferome- 
ters of Figure |§ should see ~ 100 per year. 



5. 3 Inspiral Waveforms 
They Can Bring 



and the Information 



Neutron stars and black holes have such intense self grav- 
ity that it is exceedingly difficult to deform them. Cor- 
respondingly, as they spiral inward in a compact binary, 
they do not gravitationally deform each other significantly 
until several orbits before their final coalescence |33| ^3 . 
This means that the inspiral waveforms are determined 
to high accuracy by only a few, clean parameters: the 
masses and spin angular momenta of the bodies, and the 
initial orbital elements (i.e. the elements when the waves 
enter the LIGO/VIRGO band). 

Though tidal deformations are negligible during inspi- 
ral, relativistic effects can be very important. If, for the 
moment, we ignore the relativistic effects — i.e., if we ap- 
proximate gravity as Newtonian and the wave generation 
as due to the binary's oscillating quadrupole moment ||, 
then the shapes of the inspiral waveforms h+(t) and h x (t) 
are as shown in Figure [l^. 



The left-hand graph in Figure |10| shows the waveform 
increasing in amplitude and sweeping upward in frequency 
(i.e., undergoing a "chirp") as the binary's bodies spi- 
ral closer and closer together. The ratio of the ampli- 
tudes of the two polarizations is determined by the in- 
clination 1 of the orbit to our line of sight (lower right 
in Fig. [l(]). The shapes of the individual waves, i.e. 
the waves' harmonic content, are determined by the or- 
bital eccentricity (upper right). (Binaries produced by 
normal stellar evolution should be highly circular due to 
past radiation reaction forces, but compact binaries that 
form by capture events, in dense star clusters that might 
reside in galactic nuclei |3f|, could be quite eccentric.) 
If, for simplicity, the orbit is circular, then the rate at 
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Figure 10: Waveforms from the inspiral of a compact binary, com- 
puted using Newtonian gravity for the orbital evolution and the 
quadrupole-moment approximation for the wave generation. (From 
Ref. 0.) 



which the frequency sweeps or "chirps", df/dt [or equiva- 
lently the number of cycles spent near a given frequency 
n = f 2 (df / dt)^ 1 ] is determined solely in the Newto- 
nian/quadrupole approximation, by the binary's so-called 
chirp mass, M c = (M 1 M 2 ) 3/5 /(Mi + M 2 ) 1/5 (where Mi 
and M 2 are the two bodies' masses). The amplitudes 
of the two waveforms are determined by the chirp mass, 
the distance to the source, and the orbital inclination. 
Thus (in the Newtonian/quadrupole approximation), by 
measuring the two amplitudes, the frequency sweep, and 
the harmonic content of the inspiral waves, one can de- 
termine as direct, resulting observables, the source's dis- 
tance, chirp mass, inclination, and eccentricity J3(| [57J . 

As in binary pulsar observations ||, so also here, rela- 
tivistic effects add further information: they influence the 
rate of frequency sweep and produce waveform modula- 
tions in ways that depend on the binary's dimensionless 
ratio i] = fi/M of reduced mass u = MiM 2 /(Mi + M 2 ) 
to total mass M = M\ + M 2 |40) and on the spins of 
the binary's two bodies Hj . These relativistic effects are 
reviewed and discussed at length in Refs. [^8[ Two 
deserve special mention: (i) As the waves emerge from the 
binary, some of them get backscattcrcd one or more times 
off the binary's spacetime curvature, producing wave tails. 
These tails act back on the binary, modifying its inspiral 
rate in a measurable way. (ii) If the orbital plane is in- 
clined to one or both of the binary's spins, then the spins 
drag inertial frames in the binary's vicinity (the "Lense- 
Thirring effect"), this frame dragging causes the orbit 
to precess, and the precession modulates the waveforms 
p8| , [§9[ Q . Figure [ll] shows the resulting modulation for 
a 1M Q NS spiraling into a rapidly spinning, lOM© BH. 

Remarkably, the relativistic corrections to the fre- 
quency sweep will be measurable with very high accuracy, 
even though they are typically ^ 10 per cent of the New- 
tonian contribution, and even though the typical signal 
to noise ratio will be only ~ 9 even after optimal signal 
processing. The reason is as follows E3, Ej| |3£| : 
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Figure 11: Modulational envelope for the waveform from a 1M© 
nonspinning NS spiraling into a IOMq , rapidly spinning Kerr black 
hole (spin parameter a = 1). The orbital angular momentum L 
is inclined by a = 11.3 degrees to the hole's spin angular momen- 
tum S, and the two precess around J = L + S, whose direction 
remains fixed in space as L = |L| shrinks and S = |S| = M-g^a re- 
mains constant. The precession modulates the waves by an amount 
that depends on (i) the direction to Earth (here along the initial 
L X S, i.e. out of the paper) and (ii) the orientation of the detec- 
tor's arms (here parallel to the figure's initial L and to Lx (direction 
to Earth) for and rotated 45 degrees for h x ). The figure 

shows the waveforms' modulational envelopes (in arbitrary units, 
the same for h-y. and h x ) , parametrized by the wave frequency 
/ and the number of cycles of oscillation between the indicated 
f's. The total number of precessions from / to coalescence is 
iVprec — (5/647r)(A/a//i)(7rA/q- 2 / 3 ~ 20(//10Hz)" 2 / 3 . (From 

§§[§■) 



The frequency sweep will be monitored by the method 
of "matched filters" ; in other words, the incoming, noisy 
signal will be cross correlated with theoretical templates. 
If the signal and the templates gradually get out of phase 
with each other by more than ~ 1/10 cycle as the waves 
sweep through the LIGO/VIRGO band, their cross corre- 
lation will be significantly reduced. Since the total num- 
ber of cycles spent in the LIGO/VIRGO band will be 
- 16,000 for a NS/NS binary, - 3500 for NS/BH, and 
~ 600 for BH/BH, this means that LIGO/VIRGO should 
be able to measure the frequency sweep to a fractional pre- 
cision <i 10~ 4 , compared to which the relativistic effects 
are very large. (This is essentially the same method as 
Joseph Taylor and colleagues use for high-accuracy radio- 
wave measurements of relativistic effects in binary pulsars 

§•) 

Preliminary analyses, using the theory of optimal sig- 
nal processing, predict the following typical accuracies 
for LIGO/VIRGO measurements based solely on the fre- 
quency sweep (i.e., ignoring modulational information) 
@ @ ||, 0, |§: (i) The chirp mass M c will typically 
be measured, from the Newtonian part of the frequency 
sweep, to ~ 0.04% for a NS/NS binary and - 0.3% for 
a system containing at least one BH. (ii) If we are con- 
fident (e.g., on a statistical basis from measurements of 
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many previous binaries) that the spins are a few percent 
or less of the maximum physically allowed, then the re- 
duced mass /j, will be measured to ~ 1% for NS/NS and 
NS/BH binaries, and ~ 3% for BH/BH binaries. (Here 
and below NS means a ~ 1.4M© neutron star and BH 
means a ~ 1OM black hole.) (iii) Because the frequency 
dependences of the (relativistic) \x effects and spin effects 
are not sufficiently different to give a clean separation be- 
tween fi and the spins, if we have no prior knowledge of 
the spins, then the spin//! correlation will worsen the typ- 
ical accuracy of jU by a large factor, to ~ 30% for NS/NS, 
~ 50% for NS/BH, and a factor - 2 for BH/BH g|, £wj|. 
These worsened accuracies might be improved somewhat 
by waveform modulations caused by the spin-induced pre- 
cession of the orbit |4j|, and even without modula- 
tional information, a certain combination of /j, and the 
spins will be determined to a few per cent. Much addi- 
tional theoretical work is needed to firm up the measure- 
ment accuracies. 

To take full advantage of all the information in the in- 
spiral waveforms will require theoretical templates that 
are accurate, for given masses and spins, to a fraction of a 
cycle during the entire sweep through the LIGO/VIRGO 
band. Such templates are being computed by an inter- 
national consortium of relativity theorists (Blanchet and 
Damour in France, Iyer in India, Will and Wiseman in 
the U.S., and others) p3, H, using post-Newtonian ex- 
pansions of the Einstein field equations. This enterprise 
is rather like computing the Lamb shift to high order in 
powers of the fine structure constant, for comparison with 
experiment. The terms of leading order in the mass ra- 
tio rj = n/M are being checked by a Japanese- American 
consortium (Nakamura, Sasaki, Tagoshi, Tanaka, Pois- 
son) using the Teukolsky formalism for weak perturba- 
tions of black holes |l9|, ^0|. These small- 77 calculations 
have been carried to very high post-Newtonian order for 
circular orbits and no spin s pi] , |52|, and from those re- 
sults Cutler and Flanagan |53| have estimated the order 
to which the full, finite-?? computations must be carried 
in order that systematic errors in the theoretical tem- 
plates will not significantly impact the information ex- 
tracted from the LIGO/VIRGO observational data. The 
answer appears daunting: radiation-reaction effects must 
be computed to three full post-Newtonian orders [six or- 
ders in v/c =(orbital velocity)/ (speed of light)] beyond 
the leading-order radiation reaction, which itself is 5 or- 
ders in v/c beyond the Newtonian theory of gravity. 

It is only about ten years since controversies over the 
leading-order radiation reaction [ |54| were resolved by a 
combination of theoretical techniques and binary pulsar 
observations. Nobody dreamed then that LIGO/VIRGO 
observations will require pushing post- Newtonian compu- 
tations onward from 0[(v/c) 5 ] to 0[(v/c) n ]. This require- 
ment epitomizes a major change in the field of relativity 



research: At last, 80 years after Einstein formulated gen- 
eral relativity, experiment has become a major driver for 
theoretical analyses. 

Remarkably, the goal of 0[(w/c) 11 ] is achievable. The 
most difficult part of the computation, the radiation re- 
action, has been evaluated to 0[(v/c) 9 ] beyond Newton 
by the French/Indian/ American consortium p{| and as 
of this writing, rumors have it that 0[(v / c) 10 ] is coming 
under control. 

These high-accuracy waveforms are needed only for ex- 
tracting information from the inspiral waves, after the 
waves have been discovered; they are not needed for the 
discovery itself. The discovery is best achieved using a 
different family of theoretical waveform templates, one 
that covers the space of potential waveforms in a man- 
ner that minimizes computation time instead of a man- 
ner that ties quantitatively into general relativity theory 
[38|. Such templates are in the early stage of development 

Eli m, 

LIGO/VIRGO observations of compact binary inspiral 
have the potential to bring us far more information than 
just binary masses and spins: 

• They can be used for high-precision tests of general 
relativity. In scalar-tensor theories (some of which 
are attractive alternatives to general relativity f38|), 
radiation reaction due to emission of scalar waves 
places a unique signature on the gravitational waves 
that LIGO/VIRGO would detect — a signature that 
can be searched for with high precision ]59j . 

• They can be used to measure the Universe's Hub- 
ble constant, deceleration parameter, and cosmolog- 
ical constant |3(| |37], |60[ [il) . The keys to such mea- 
surements are that: (i) Advanced interferometers in 
LIGO/VIRGO will be able to see NS/NS out to cos- 
mological redshifts z ~ 0.3, and NS/BH out to z ~ 2. 
(ii) The direct observables that can be extracted from 
the observed waves include the source's luminosity 
distance ?*l (measured to accuracy ~ 10 per cent in 
a large fraction of cases), and its direction on the sky 
(to accuracy ~ 1 square degree) — accuracies good 
enough that only one or a few electromagnetically- 
observed clusters of galaxies should fall within the 3- 
dimensional gravitational error boxes, thereby giving 
promise to joint gravitational/electromagnetic statis- 
tical studies, (iii) Another direct gravitational ob- 
servable is (1 + z)M where z is redshift and M is 
any mass in the system (measured to the accuracies 
quoted above). Since the masses of NS's in bina- 
ries seem to cluster around IAM@, measurements of 
(1 + z)M can provide a handle on the redshift, even 
in the absence of electromagnetic aid. 

• For a NS or small BH spiraling into a massive ~ 50 to 
500M Q BH, the inspiral waves will carry a "map" of 
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the spacetime geometry around the big hole — a map 
that can be used, e.g., to test the the orem that "a 
black hole has no hair" Ip2fl ; cf . Section |8.3| below. 



5.4 Coalescence Waveforms and their Informa- 
tion 

The waves from the binary's final coalescence can bring 
us new types of information. 

BH/BH Coalescence 

In the case of a BH/BH binary, the coalescence will ex- 
cite large-amplitude, highly nonlinear vibrations of space- 
time curvature near the coalescing black-hole horizons — a 
phenomenon of which we have very little theoretical un- 
derstanding today. Especially fascinating will be the case 
of two spinning black holes whose spins are not aligned 
with each other or with the orbital angular momentum. 
Each of the three angular momentum vectors (two spins, 
one orbital) will drag space in its vicinity into a tornado- 
like swirling motion — the general relativistic "dragging of 
inertial frames," so the binary is rather like two torna- 
dos with orientations skewed to each other, embedded in- 
side a third, larger tornado with a third orientation. The 
dynamical evolution of such a complex configuration of 
coalescing spacetime warpage (as revealed by its emitted 
waves) might bring us surprising new insights into rel- 
ativistic gravity [[14j. Moreover, if the sum of the BH 
masses is fairly large, ~ 40 to 200M©, then the waves 
should come off in a frequency range / ~ 40 to 200 
Hz where the LIGO/VIRGO broad-band interferometers 
have their best sensitivity and can best extract the infor- 
mation the waves carry. 

To get full value out of such wave observations will re- 
quire f(33[l having theoretical computations with which to 
compare them. There is no hope to perform such compu- 
tations analytically; they can only be done as supercom- 
puter simulations. The development of such simulations 
is being pursued by several research groups, including an 
eight-university American consortium of numerical rela- 
tivists and computer scientists called the Two-Black-Holc 
Grand Challenge Alliance [Q (Co-PIs: Richard Matzner 
and Jim Browne, U. Texas Austin; Larry Smarr, Ed Sei- 
del, Paul Saylor, Faisal Saied, U. Illinois Urbana; Geoffrey 
Fox, Syracuse U.; Stu Shapiro and Saul Teukolsky, Cornell 
U.; Jim York and Charles Evans, U. North Carolina; Sam 
Finn, Northwestern U.; Pablo Laguna, Pennsylvania State 
U.; and Jeff Winicour, U. Pittsburgh). I have a bet with 
Matzner, the lead PI of this alliance, that LIGO/VIRGO 
will discover waves from such coalescences with misaligned 
spins before the Alliance is able to compute them. 

NS/NS Coalescence 



The final coalescence of NS /NS binaries should produce 
waves that are sensitive to the equation of state of nu- 
clear matter, so such coalescences have the potential to 
teach us about the nuclear equation of state [jl4| |38|. In 
essence, LIGO/VIRGO will be studying nuclear physics 
via the collisions of atomic nuclei that have nucleon num- 
bers A ~ 10 57 — somewhat larger than physicists are nor- 
mally accustomed to. The accelerator used to drive these 
nuclei up to the speed of light is the binary's self gravity, 
and the radiation by which the details of the collisions are 
probed is gravitational. 

Unfortunately the final NS/NS coalescence will emit its 
gravitational waves in the kHz frequency band (800Hz <, 
f & 2500Hz) where photon shot noise will prevent them 
from being studied by the standard, "workhorse," broad- 
band interferometers of Figure |^. However, a specially 
configured ("dual-recycled") interferometer invented by 
Brian Meers |65| , which could have enhanced sensitivity 
in the kHz region at the price of reduced sensitivity else- 
where, may be able to measure the waves and extract 
their equation of state information, as might massive, 
spherical, resonant-mass detectors 38, ^36|. Such mea- 



surements will be very difficult and are likely only when 
the LIGO/VIRGO network has reached a mature stage 
A number of research groups (64 pq, j 



are 

engaged in numerical astrophysics simulations of NS/NS 
coalescence, with the goal not only to predict the emitted 
gravitational waveforms and their dependence on equation 
of state, but also (more immediately) to learn whether 
such coalescences might power the 7-ray bursts that have 
been a major astronomical puzzle since their discovery in 
the early 1970s. 

NS/NS coalescence is currently a popular explanation 
for the 7-ray bursts because (i) the bursts are isotropi- 
cally distributed on the sky, (ii) they have a distribution 
of number versus intensity that suggests they might lie 
at near-cosmological distances, and (iii) their event rate 
is roughly the same as that predicted for NS/NS coales- 
cence (~ 1000 per year out to cosmological distances, if 
they are cosmological). If LIGO/VIRGO were now in 
operation and observing NS/NS inspiral, it could report 
definitively whether or not the 7-bursts are produced by 
NS/NS binaries; and if the answer were yes, then the 
combination of 7-burst data and gravitational-wave data 
could bring valuable information that neither could bring 
by itself. For example, it would reveal when, to within a 
few msec, the 7-burst is emitted relative to the moment 
the NS's first begin to touch; and by comparing the 7 
and gravitational times of arrival, we might test whether 
gravitational waves propagate with the speed of light to 
a fractional precision of ~ 0.01sec/3 x 10 9 lyr = 10 -19 . 

NS/BH Coalescence 
A NS spiraling into a BH of mass M <; 10M Q should 
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be swallowed more or less whole. However, if the BH is 
less massive than roughly 10M©, and especially if it is 
rapidly rotating, then the NS will tidally disrupt before 
being swallowed. Little is known about the disruption 
and accompanying waveforms. To model them with any 
reliability will likely require full numerical relativity, since 
the circumferences of the BH and NS will be comparable 
and their physical separation at the moment of disruption 
will be of order their separation. As with NS/NS, the coa- 
lescence waves should carry equation of state information 
and will come out in the kHz band, where their detection 
will require advanced, specialty detectors. 

Christodoulou Memory 

As the coalescence waves depart from their source, their 
energy creates (via the nonlinearity of Einstein's field 
equations) a secondary wave called the "Christodoulou 
memory j7|, fr| ||. Whereas the primary waves may 
have frequencies in the kHz band, the memory builds 
up on the timescale of the primary energy emission pro- 
file, which is likely to be of order 0.01 sec, correspond- 
ing to a memory frequency in the optimal band for the 
LIGO/VIRGO workhorse interferometers, - 100Hz. Un- 
fortunately, the memory is so weak that only very ad- 
vanced interferometers have much chance of detecting and 
studying it — and then, perhaps only for BH/BH coales- 
cences and not for NS/NS or NS/BH j7§. 

§6 Other High-Prequency Sources 

6.1 Stellar Core Collapse and Supernovae 

When the core of a massive star has exhausted its sup- 
ply of nuclear fuel, it collapses to form a neutron star or 
black hole. In some cases, the collapse triggers and pow- 
ers a subsequent explosion of the star's mantle — a super- 
nova explosion. Despite extensive theoretical efforts for 
more than 30 years, and despite wonderful observational 
data from Supernova 1987A, theorists are still far from a 
definitive understanding of the details of the collapse and 
explosion. The details are highly complex and may differ 
greatly from one core collapse to another J7(J . 

Several features of the collapse and the core's subse- 
quent evolution can produce significant gravitational ra- 
diation in the high-frequency band. We shall consider 
these features in turn, the most weakly radiating first. 

Boiling of the Newborn Neutron Star 

Even if the collapse is spherical, so it cannot radiate 
any gravitational waves at all, it should produce a con- 
vectively unstable neutron star that "boils" vigorously 
(and nonsphcrically) for the first ~ 0.1 second of its 
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Figure 12: Characteristic amplitudes of the gravitational waves 
from various processes accompanying stellar core collapse and su- 
pernovae, compared with projected sensitivities of LIGO's interfer- 
ometers. 



life |77[ |. The boiling dredges up high-temperature nu- 
clear matter (T ~ 10 12 K) from the neutron star's cen- 
tral regions, bringing it to the surface (to the "neutrino- 
sphere"), where it cools by neutrino emission before be- 
ing swept back downward and reheated. Burrows esti- 
mates 1 78, [79| that the boiling should generate n ~ 10 
cycles of gravitational waves with frequency / ~ 100Hz 
and amplitude h ~ 3 x 10~ 22 (30kpc/r) (where r is the 



distance to the source), corresponding to a characteristic 
amplitude h c ~ h^/n - 10~ 21 (30kpc/r); c f. Figure @ 
LIGO/VIRGO will be able to detect such waves only in 
the local group of galaxies, where the supernova rate is 
probably no larger than ~ 1 each 10 years. However, neu- 
trino detectors have a similar range, and there could be a 
high scientific payoff from correlated observations of the 
gravitational waves emitted by the boiling's mass motions 
and neutrinos emitted from the boiling neutrino-sphere. 

Axisymmetric Collapse, Bounce, and Oscillations 

Rotation will centrifugally flatten the collapsing core, 
enabling it to radiate as it implodes. If the core's angular 
momentum is small enough that centrifugal forces do not 
halt or strongly slow the collapse before it reaches nuclear 
densities, then the core's collapse, bounce, and subsequent 
oscillations are likely to be axially symmetric. Numerical 
simulations |30| [y] show that in this case the waves from 
collapse, bounce, and oscillation will be quite weak: the 
total energy radiated as gravitational waves is not likely 
to exceed ~ 10~ 7 solar masses (about 1 part in a million 
of the collapse energy) and might often be much less than 
this; and correspondingly, the waves' characteristic am- 
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plitude will be h c <S 3 x 10 _21 (30kpc/r). These collapse- 
and-bounce waves will come off at frequencies ~ 200 Hz 
to ~ 1000 Hz, and will precede the boiling waves by a 
fraction of a second. Like the boiling waves, they prob- 
ably cannot be seen by LIGO/ VIRGO beyond the local 
group of galaxies and thus will be a very rare occurrence. 

Rotation- Induced Bars and Break- Up 

If the core's rotation is large enough to strongly flat- 
ten the core before or as it reaches nuclear density, then 
a dynamical and/or secular instability is likely to break 
the core's axisymmetry. The core will be transformed 
into a bar-like configuration that spins end-over-end like 
an American football, and that might even break up into 
two or more massive pieces. In this case, the radiation 
from the spinning bar or orbiting pieces could be almost 
as strong as that from a coalescing neutron-star binary, 
and thus could be seen by the LIGO/ VIRGO first inter- 
ferometers out to the distance of the Virgo cluster (where 
the supernova rate is several per year) and by advanced 
interferometers out to several hundred Mpc (supernova 
rate ~ 10 4 per year); cf. Figure [jj]. It is far from clear 
what fraction of collapsing cores will have enough angular 
momentum to break their axisymmetry, and what fraction 
of those will actually radiate at this high rate; but even if 
only ~ 1/1000 or 1/10 4 do so, this could ultimately be a 
very interesting source for LIGO/ VIRGO. 

Several specific scenarios for such non-axisymmetry 
have been identified: 

Centrifugal hangup at ~ 100km radius: If the pre- 
collapse core is rapidly spinning (e.g., if it is a white dwarf 
that has been spun up by accretion from a companion), 
then the collapse may produce a highly flattened, cen- 
trifugally supported disk with most of its mass at radii 
R ~ 100km, which then (via instability) may transform 
itself into a bar or may bifurcate. The bar or bifurcated 
lumps will radiate gravitational waves at twice their ro- 
tation frequency, / ~ 100Hz — the optimal frequency 
for LIGO/VIRGO interferometers. To shrink on down to 
~ 10km size, this configuration must shed most of its an- 
gular momentum. If & substantial fraction of the angular 
momentum goes into gravitational waves, then indepen- 
dently of the strength of the bar, the waves will be nearly 
as strong as those from a coalescing binary. The rea- 
son is this: The waves' amplitude h is proportional to the 
bar's ellipticity e, the number of cycles n of wave emission 
is proportional to 1/e 2 , and the characteristic amplitude 
h c = hsfn is thus independent of the ellipticity and is 
about the same whether the configuration is a bar or is 
two lumps The resulting waves will thus have h c 

roughly half as large, at / ~ 100Hz, as the h c from a 
NS /NS binary (half as large because each lump might be 
half as massive as a NS), and the waves will chirp upward 
in frequency in a manner similar to those from a binary. 



It is rather likely, however, that most of excess angu- 
lar momentum does not go into gravitational waves, but 
instead goes largely into hydrodynamic waves as the bar 
or lumps, acting like a propeller, stir up the surrounding 
stellar mantle. In this case, the radiation will be corre- 
spondingly weaker. 

Centrifugal hangup at ~ 20km radius: Lai and 
Shapiro J82fl have explored the case of centrifugal hangup 
at radii not much larger than the final neutron star, say 
R ~ 20km. Using compressible ellipsoidal models, they 
have deduced that, after a brief period of dynamical bar- 
mode instability with wave emission at / ~ 1000Hz (ex- 
plored by Houser, Centrella, and Smith |83]|), the star 
switches to a secular instability in which the bar's angu- 
lar velocity gradually slows while the material of which 
it is made retains its high rotation speed and circu- 
lates through the slowing bar. The slowing bar emits 
waves that sweep downward in frequency through the 
LIGO/VIRGO optimal band / ~ 100Hz, toward ~ 10Hz. 
The characteristic amplitude (Fig. [ljf) is only modestly 
smaller than for the upward-sweeping waves from hangup 
at R ~ 100km, and thus such waves should be detectable 
near the Virgo Cluster by the first LIGO/VIRGO inter- 
ferometers, and at distances of a few lOOMpc by advanced 
interferometers. 

Successive fragmentations of an accreting, new- 
born neutron star: Bonnell and Pringle ||3] have fo- 
cused on the evolution of the rapidly spinning, newborn 
neutron star as it quickly accretes more and more mass 
from the pre-supernova star's inner mantle. If the accret- 
ing material carries high angular momentum, it may trig- 
ger a renewed bar formation, lump formation, wave emis- 
sion, and coalescence, followed by more accretion, bar and 
lump formation, wave emission, and coalescence. Bon- 
nell and Pringle speculate that hydrodynamics, not wave 
emission, will drive this evolution, but that the total en- 
ergy going into gravitational waves might be as large as 
~ 10 _3 A/ o . This corresponds to h c ~ 10~ 21 (10Mpc/r). 

6.2 Spinning Neutron Stars; Pulsars 

As the neutron star settles down into its final state, its 
crust begins to solidify (crystalize) . The solid crust will 
assume nearly the oblate axisymmetric shape that cen- 
trifugal forces are trying to maintain, with poloidal ellip- 
ticity e p oc(angular velocity of rotation) 2 . However, the 
principal axis of the star's moment of inertia tensor may 
deviate from its spin axis by some small "wobble angle" 
9 W , and the star may deviate slightly from axisymmetry 
about its principal axis; i.e., it may have a slight ellipticity 
e e -C e p in its equatorial plane. 

As this slightly imperfect crust spins, it will radiate 
gravitational waves |8f| : e e radiates at twice the rotation 
frequency, / = 2/ rot with h cx e e , and the wobble angle 



couples to e p to produce waves at / = / ro t + / pr cc (the 
precessional sideband of the rotation frequency) with am- 
plitude h oc 9 w e p . For typical neutron-star masses and 
moments of inertia, the wave amplitudes are 

h „ 6 x IP' 25 ( J22-) 2 ( I?**) h ° r 9 f P ) ■ (5) 
V 500Hz / \ r J \ 10- 6 J V 1 

The neutron star gradually spins down, due in part 
to gravitational-wave emission but perhaps more strongly 
due to electromagnetic torques associated with its spin- 
ning magnetic field and pulsar emission. This spin-down 
reduces the strength of centrifugal forces, and thereby 
causes the star's poloidal ellipticity e p to decrease, with an 
accompanying breakage and resolidification of its crust's 
crystal structure (a "starquake") |8^| . In each starquake, 
9 W , e e , and e p will all change suddenly, thereby changing 
the amplitudes and frequencies of the star's two gravita- 
tional "spectral lines" / = 2/ rot and / = /rot + /prcc- After 
each quake, there should be a healing period in which the 
star's fluid core and solid crust, now rotating at differ- 
ent speeds, gradually regain synchronism. By monitor- 
ing the amplitudes, frequencies, and phases of the two 
gravitational- wave spectral lines, and by comparing with 
timing of the electromagnetic pulsar emission, one might 
learn much about the physics of the neutron-star interior. 

How large will the quantities e e and 9 w e p be? Rough es- 
timates of the crustal shear moduli and breaking strengths 
suggest an upper limit in the range e max ~ 1CT 4 to 10~ 6 , 
and it might be that typical values are far below this. 
We are extremely ignorant, and correspondingly there is 
much to be learned from searches for gravitational waves 
from spinning neutron stars. 

One can estimate the sensitivity of LIGO/VIRGO (or 
any other broad-band detector) to the periodic waves from 
such a source by multiplying the waves' amplitude h by 
the square root of the number of cycles over which one 
might integrate to find the signal, n — ff where f is the 
integration time. The resulting effective signal strength, 
hy/n, is larger than h by 

= \fj? = 10 5 ( — - — 1 ' ( 1 1 ■ (6) 

This h^/n should be compared (i) to the detector's rms 
broad-band noise level for sources in a random direction, 
Vbh lms , to deduce a signal-to-noise ratio, or (ii) to Iisb to 
deduce a sensitivity for high-confidence detection when 
one does not know the waves' frequency in advance 
Such a comparison suggests that the first interferometers 
in LIGO/VIRGO might possibly see waves from nearby 
spinning neutron stars, but the odds of success are very 
unclear. 

The deepest searches for these nearly periodic waves 
will be performed by narrow-band detectors, whose sen- 



sitivities are enhanced near some chosen frequency at the 
price of sensitivity loss elsewhere — e.g., dual-recycled in- 
terferometers |H or resonant-mass antennas (Section |]) . 
With "advanced-detector technology," dual-recycled in- 
terferometers might be able to detect with confidence all 
spinning neutron stars that have j(| 

(E . or9 . E ^ 3xl0 -.«(^) 2 ( T5 i 5 _) 2 . P ) 

There may well be a large number of such neutron stars in 
our galaxy; but it is also conceivable that there are none. 
We are extremely ignorant. 

Some cause for optimism arises from several physi- 
cal mechanisms that might generate radiating ellipticities 
large compared to 3 x 10~ 10 : 

• It may be that, inside the superconducting cores of 
many neutron stars, there are trapped magnetic fields 
with mean strength B cmc ~ 10 13 G or even 10 15 G. 
Because such a field is actually concentrated in flux 
tubes with B = B CT a ~ 6 x 10 14 G surrounded by 
field-free superconductor, its mean pressure is ps = 
B cole B CI it/8Tr. This pressure could produce a radiat- 
ing ellipticity e c ~ 9 w e p ~ p B /p ~ 10" 8 B CO r C /10 13 G 
(where p is the core's material pressure). 

• Accretion onto a spinning neutron star can drive 
precession (keeping 9 W substantially nonzero), and 
thereby might produce measurably strong waves | g7J . 

• If a neutron star is born rotating very rapidly, then 
it may experience a gravitational-radiation-reaction- 
driven instability. In this "CFS" (Chandrasekhar, 
[p8| Friedman, Schutz |||) instability, density waves 
propagate around the star in the opposite direction to 
its rotation, but are dragged forward by the rotation. 
These density waves produce gravitational waves that 
carry positive energy as seen by observers far from the 
star, but negative energy from the star's viewpoint; 
and because the star thinks it is losing negative en- 
ergy, its density waves get amplified. This intriguing 
mechanism is similar to that by which spiral density 
waves are produced in galaxies. Although the CFS 
instability was once thought ubiquitous for spinning 
stars p9| , pC[ |, we now know that neutron-star vis- 
cosity will kill it, stabilizing the star and turning off 
the waves, when the star's temperature is above some 
limit - 10 10 K |9lJ and below some limit - 10 9 K @; 
and correspondingly, the instability should operate 
only during the first few years of a neutron star's life, 
when 10 9 K £ T £ 10 10 K. 



6.3 Stochastic Background 

There should be a stochastic background of gravitational 
waves in the high-frequency band produced by processes 
in the early universe. Because this background will ex- 
tend over all gravitational-wave frequencies, not just high 
frequencies, we shall delay discussing it until Section ^|. 

§7 LISA: The Laser Interferometer 
Space Antenna 

Turn, now, from the high-frequency band, 1 — 10 Hz, 
to the low-frequency band, 10~ 4 — 1 Hz. At present, 
the most sensitive gravitational-wave searches at low fre- 
quencies are those carried out by researchers at NASA's 
Jet Propulsion Laboratory, using microwave-frequency 
Doppler tracking of interplanetary spacecraft. These 
searches are done at rather low cost, piggy-back on mis- 
sions designed for other purposes. Although they have a 
possibility of success, the odds are against them. Their 
best past sensitivities to bursts, for example, have been 
/isb ~ 10~ 14 , and prospects are good for reaching ~ 
10 -15 — 10 -16 in the next 5 to 10 years. However, the 
strongest low-frequency bursts arriving several times per 
year might be no larger than ~ 10~ 18 ; and the domain of 
an assured plethora of signals is Iisb ~ 10~ 19 — 10~ 20 . 

To reach into this assured-detection domain will almost 
certainly require switching from microwave-frequency 
tracking of spacecraft (with its large noise due to fluc- 
tuating dispersion in the troposphere and interplanetary 
plasma) to optical tracking. Such a switch is planned for 
the 2014 time frame or sooner, when the European Space 
Agency (ESA) and/or NASA is likely to fly the Laser In- 
terferometer Space Antenna, LISA. 

7.1 Mission Status 

LISA is largely an outgrowth of 15 years of studies by 
Peter Bender and colleagues at the University of Col- 
orado. In 1990, NASA's Ad Hoc Committee on Grav- 
itation Physics and Astronomy selected a LISA-type 
gravitational-wave detector as its highest priority in the 
large space mission category |93j ; and since then enthusi- 
asm for LISA has continued to grow within the American 
gravitation community. Unfortunately, the prospects for 
NASA to fly such a mission did not look good in the early 
1990s. By contrast, prospects in Europe looked much bet- 
ter, so a largely European consortium was put together in 
1993, under the leadership of Karsten Danzmann (Han- 
nover) and James Hough (Glasgow), to propose LISA 
to the European Space Agency. This proposal has met 
with considerable success; LISA might well achieve ap- 
proval to fly as an ESA Cornerstone Mission around 2014 
p4|. Members of the American gravitation community 




Figure 13: LISA'S orbital configuration, with LISA magnified in arm 
length by a factor ~ 10 relative to the solar system. 

and members of the LISA team hope that NASA will join 
together with ESA in this endeavor, and that working 
jointly, ESA and NASA will be able to fly LISA consider- 
ably sooner than 2014. 

7.2 Mission Configuration 

As presently conceived, LISA will consist of six compact, 
drag-free spacecraft (i.e. spacecraft that are shielded from 
buffeting by solar wind and radiation pressure, and that 
thus move very nearly on geodesies of spacetime). All 
six spacecraft would be launched simultaneously by a sin- 
gle Ariane rocket. They would be placed into the same 
heliocentric orbit as the Earth occupies, but would fol- 
low 20° behind the Earth; cf. Figure [l^. The spacecraft 
would fly in pairs, with each pair at the vertex of an equi- 
lateral triangle that is inclined at an angle of 60° to the 
Earth's orbital plane. The triangle's arm length would 
be 5 million km (10 6 times larger than LIGO's arms!). 
The six spacecraft would track each other optically, us- 
ing one- Watt YAG laser beams. Because of diffraction 
losses over the 5 x 10 6 km arm length, it is not feasible 
to reflect the beams back and forth between mirrors as 
is done with LIGO. Instead, each spacecraft will have its 
own laser; and the lasers will be phase locked to each 
other, thereby achieving the same kind of phase-coherent 
out-and-back light travel as LIGO achieves with mirrors. 
The six-laser, six-spacecraft configuration thereby func- 
tions as three, partially independent but partially redun- 
dant, gravitational-wave interferometers. 

7.3 Noise and Sensitivity 

Figure |l4| depicts the expected noise and sensitivity of 
LISA in the same language as we have used for LIGO 
(Fig. . The curve at the bottom of the stippled region is 
forms > the rms noise, in a bandwidth equal to frequency, for 
waves with optimum direction and polarization. The top 
of the stippled region is hsB — 5v5/irms> the sensitivity 
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Figure 14: LISA'S projected broad-band noise /i rms and sensitivity 
to bursts /iSBi compared with the strengths of the waves from several 
low-frequency sources. [Note: When members of the LISA team 
plot curves analogous to this, they show the sensitivity curve (top 
of stippled region) in units of the amplitude of a periodic signal 
that can be detected with S/N = 5 in one year of integration; that 
sensitivity to periodic sources is related to the hgB used here by 

h S p = hss/y/f ■ 3 x 10W] 

for high-confidence detection (S/N = 5) of a broad-band 
burst coming from a random direction, assuming Gaussian 
noise. 

At frequencies / ^ 10~ 3 Hz, LISA's noise is due to pho- 
ton counting statistics (shot noise) . The noise curve steep- 
ens at / ~ 3 x 10~ 2 Hz because at larger / than that, the 
waves' period is shorter than the round-trip light travel 
time in one of LISA's arms. Below 10 _3 Hz, the noise is 
due to buffeting-induced random motions of the space- 
craft that are not being properly removed by the drag- 
compensation system. Notice that, in terms of dimen- 
sionless amplitude, LISA's sensitivity is roughly the same 
as that of LIGO's first interferometers (Fig. |^), but at 
100,000 times lower frequency. Since the waves' energy 
flux scales as f 2 h 2 , this corresponds to 10 10 better energy 
sensitivity than LIGO. 



7.4 Observational Strategy 

LISA can detect and study, simultaneously, a wide vari- 
ety of different sources scattered over all directions on the 
sky. The key to distinguishing the different sources is the 
different time evolution of their waveforms. The key to 
determining each source's direction, and confirming that 
it is real and not just noise, is the manner in which its 
waves' amplitude and frequency are modulated by LISA's 



complicated orbital motion — a motion in which the inter- 
ferometer triangle rotates around its center once per year, 
and the interferometer plane precesses around the normal 
to the Earth's orbit once per year. Most sources will be 
observed for a year or longer, thereby making full use of 
these modulations. 

§8 Low-Frequency Gravitational- Wave 
Sources 

8. 1 Waves from Binary Stars 

LISA has a large class of guaranteed sources: short-period 
binary stars in our own galaxy. A specific example is the 
classic binary 44 i Boo (HD133640), a 1.35M /O.68M 
system just 12 parsecs from Earth, whose wave frequency 
/ and characteristic amplitude h c = h^/n are depicted in 
Figure [l4|. (Here h is the waves' actual amplitude and 
n = ff is the number of wave cycles during f =1 year of 
signal integration). Since 44 i Boo lies right on the hsB 
curve, its signal to noise ratio in one year of integration 
should be S/N = 5. 

To have an especially short period, a binary must be 
made of especially compact bodies — white dwarfs (WD), 
neutron stars (NS), and/or black holes (BH). WD/WD 
binaries are thought to be so numerous that they may 
produce a stochastic background of gravitational waves, 
at the level shown in Figure O, that will hide some other 
interesting waves from view Ba j. Since WD/ WD binaries 
are very dim optically, their actual numbers are not known 
for sure; Figure ^ might be an overestimate. 

Assuming a NS/NS coalescence rate of 1 each 10 5 years 
in our galaxy |t], ||, the shortest period NS/NS bi- 
nary should have a remaining life of about 5 x 10 years, 
corresponding to a gravitational-wave frequency today of 
/ ~ 5 x 10 _3 Hz, an amplitude (at about lOkpc distance) 
ft~4x 10 -22 , and a characteristic amplitude (with one 
year of integration time) h c ~ 2 x 10 -19 . This is de- 
picted in Figure [l4] at the right edge of the region marked 
"brightest NS/NS binaries". These brightest NS/NS bi- 
naries can be studied by LISA with the impressive signal 
to noise ratios S/N ~ 50 to 500. 

8.2 Waves from the Coalescence of Massive 
Black Holes in Distant Galaxies 

LISA would be a powerful instrument for studying mas- 
sive black holes in distant galaxies. Figure |lj shows, as 
examples, the waves from several massive black hole bina- 
ries at 3Gpc distance from Earth (a cosmological rcdshift 
of unity). The waves sweep upward in frequency (right- 
ward in the diagram) as the holes spiral together. The 
black dots show the waves' frequency one year before the 



holes' final collision and coalescence, and the arrowed lines 
show the sweep of frequency and characteristic amplitude 
h c = h^/ri during that last year. For simplicity, the fig- 
ure is restricted to binaries with equal-mass black holes: 
1O 4 M /1O 4 M , lO 5 Af /lO 5 M , and 1O 6 M Q /1O 6 M . 

By extrapolation from these three examples, we see that 
LISA can study much of the last year of inspiral, and the 
waves from the final collision and coalescence, whenever 
the holes' masses are in the range 3 x 1O 4 M < M < 
1O 8 M . Moreover, LISA can study the final coalescences 
with remarkable signal to noise ratios: S/N > 1000. Since 
these are much larger S/N's than LIGO/ VIRGO is likely 
to achieve, we can expect LISA to refine the experimental 
understanding of black-hole physics, and of highly nonlin- 
ear vibrations of warped spacetime, which LIGO/VIRGO 
initiates — provided the rate of massive black-hole coales- 
cences is of order one per year in the Universe or higher. 
The rate might well be that high, but it also might be 
much lower. 

By extrapolating Figure [lj to lower BH/BH masses, 
we see that LISA can observe the last few years of inspi- 
ral, but not the final collisions, of binary black holes in 
the range 100Af Q <, M ^ 1O 4 M , out to cosmological 
distances. 

Extrapolating the BH/BH curves to lower frequencies 
using the formula (time to final coalescence) cx /~ 8 ^ 3 , 
we see that equal-mass BH/BH binaries enter LISA's fre- 
quency band roughly 1000 years before their final coa- 
lescences, more or less independently of their masses, for 
the range 100Af Q <, M <, 10 6 Af Q . Thus, if the coales- 
cence rate were to turn out to be one per year, LISA 
would see roughly 1000 additional massive binaries that 
are slowly spiraling inward, with inspiral rates df / dt read- 
ily measurable. From the inspiral rates, the amplitudes 
of the two polarizations, and the waves' harmonic con- 
tent, LISA can determine each such binary's luminosity 
distance, redshifted chirp mass (1 + z)M c , orbital incli- 
nation, and eccentricity; and from the waves' modulation 
by LISA's orbital motion, LISA can learn the direction to 
the binary with an accuracy of order one degree. 

8.3 Waves from Compact Bodies Spiraling into 
Massive Black Holes in Distant Galaxies 

When a compact body with mass fj, spirals into a much 
more massive black hole with mass M, the body's orbital 
energy E at fixed frequency / (and correspondingly at 
fixed orbital radius a) scales as E cx /i, the gravitational- 
wave luminosity E scales as E cx /i 2 , and the time to 
final coalescence thus scales as t ~ E/E cx l//i. This 
means that the smaller is fi/M, the more orbits are spent 
in the hole's strong-gravity region, a <^ lOGM/c 2 , and 
thus the more detailed and accurate will be the map of 



the hole's spacetime geometry, which is encoded in the 
emitted waves. 

For holes observed by LIGO/VIRGO, the most extreme 
mass ratio that we can hope for is fi/M ~ lAf /3OOAf , 
since for M > 3OOM the inspiral waves are pushed to 
frequencies below the LIGO/VIRGO band. This limit 
on fi/M seriously constrains the accuracy with which 
LIGO/VIRGO can hope to map out the spacetime ge- 
ometries of black holes and test the black-hole no-hair 
theorem |6^] (end of Section 5.3). By contrast, LISA can 
observe the final inspiral waves from objects of any mass 
M ^ O.5M spiraling into holes of mass 3 x 1O 5 M < 
M <, 3 x 1O 7 M . 

Figure |l4| shows the example of a 1OM black hole spi- 
raling into a lO 6 Af hole at 3Gpc distance. The inspi- 
ral orbit and waves are strongly influenced by the hole's 
spin. Two cases are shown fl96| : an inspiraling circular 
orbit around a non-spinning hole, and a prograde, circu- 
lar, equatorial orbit around a maximally spinning hole. 
In each case the dot at the upper left end of the arrowed 
curve is the frequency and characteristic amplitude one 
year before the final coalescence. In the nonspinning case, 
the small hole spends its last year spiraling inward from 
r ~ TAGM/c 2 (3.7 Schwarzschild radii) to its last sta- 
ble circular orbit at r — 6GM/c 2 (3 Schwarzschild radii). 
In the maximal spin case, the last year is spent travel- 
ing from r = QGM/c 2 (3 Schwarzschild radii) to the last 
stable orbit at r — GM/c 2 (half a Schwarzschild radius). 
The ~ 10 5 cycles of waves during this last year should 
carry, encoded in themselves, rather accurate values for 
the massive hole's lowest few multipole moments p7| . If 
the measured moments satisfy the "no-hair" theorem (i.e., 
if they are all determined uniquely by the measured mass 
and spin in the manner of the Kerr metric) , then we can 
be sure the central body is a black hole. If they violate 
the no-hair theorem, then (assuming general relativity is 
correct), either the central body was not a black hole, or 
an accretion disk or other material was perturbing its or- 
bit [68| . From the evolution of the waves one can hope to 
determine which is the case, and to explore the properties 
of the central body and its environment f32f . 

Models of galactic nuclei, where massive holes reside, 
suggest that inspiraling stars and small holes typically 
will be in rather eccentric orbits J99| . This is because they 
get injected into such orbits via gravitational deflections 
off other stars, and by the time gravitational radiation 
reaction becomes the dominant orbital driving force, there 
is not enough inspiral left to fully circularize their orbits. 
Such orbital eccentricity will complicate the waveforms 
and complicate the extraction of information from them. 
Efforts to understand the emitted waveforms are just now 
getting underway. 

The event rates for inspiral into massive black holes are 
not at all well understood. However, since a significant 
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fraction of all galactic nuclei are thought to contain mas- 
sive holes, and since white dwarfs and neutron stars, as 
well as small black holes, can withstand tidal disruption as 
they plunge toward the massive hole's horizon, and since 
LISA can see inspiraling bodies as small as ~ 0.5M Q out 
to 3Gpc distance, the event rate is likely to be interest- 
ingly large. 



somewhat beyond. The observational limit at 10 18 H z 
implies that this constant value is fl g <J 3x 10~ 14 [Il2, 102|. 



\9 The Stochastic 
Background 



Gravitational- Wave 



Processes in the early universe should have produced a 
stochastic background of gravitational waves that extends 
through the entire frequency range from extremely low 
frequencies / ~ 10~ 18 Hz to the high-frequency band / ~ 
1 — 10 4 Hz and beyond. 

9.1 Primordial Gravitational Waves 

The most interesting background would be that produced 
in the big bang itself. Zel'dovich and Novikov have esti- 



mated [10C] that the optical thickness of primordial mat- 
ter to gravitational waves has been small compared to 
unity at all times since the Planck era, when space and 
time came into being, and that therefore primordial grav- 
itational waves (by contrast with electromagnetic) should 
not have been thermalized by interactions with matter. 
On the other hand (as Grishchuk has shown [101]), what- 
ever might have been the state of the graviton field when 
it emerged from the big bang's Planck era, it should have 
interacted with the subsequent, early-time expansion of 
the universe to produce, via parametric amplification, a 
rich spectrum of stochastic waves today. The details of 
that spectrum depend on what emerged from the Planck 
era and on the evolution a(t) of the universal expansion 
factor at early times. 

The gravitational- wave spectrum is generally described 
by the quantity £l g (f) = (energy density in a bandwidth 
equal to frequency /)/ (energy density required to close 
the universe); cf. Section |2.3| . The observed quadrupo- 
lar anisotropy of the cosmic microwave radiation places 



a limit fl g <, 10~ 9 at / - 10~ 18 Hz (Section |2J). It is 



fashionable to extrapolate this limit to higher frequen- 
cies by assuming that the graviton field emerged from the 
Planck era in its vacuum state, and assuming that the 
universal expansion a(t) was that of an inflationary era 
a oc e Ht for some constant H, followed by a radiation- 
dominated Friedman era a oc i 1 / 2 , followed by the present 
matter-dominated era a oc i 2 / 3 . This standard model 
produces a flat spectrum Q g independent of / for all 
waves that entered our cosmological horizon during the 
radiation-dominated era, which means at all frequencies 
from ~ 10~ 16 Hz up through the high-frequency band and 



So weak a background cannot be detected by LIGO in the 
high-frequency band, nor by LISA at low frequencies, nor 
by pulsar timing at very low frequencies. LIGO's limit- 
ing sensitivities will correspond to £l g ~ (a few) x 10~ 7 
at / ~ 10 2 Hz for the first interferometers, and fl g ~ 
(a few) x 10~ 10 for advanced interferometers [103 ; LISA's 
sensitivity will correspond to £l g ~ (a few) x 10 _ 10 at 
/ ~ 10 -3 Hz; and the present pulsar timing measurements 
correspond to fl g ~ (a few) x 10~ 8 at / ~ 4 x 10~ 9 Hz 



(Section and Ref. |9|). 

On the other hand, if the graviton field did not begin 
in its vacuum state, or if the equation of state in the 
very early Friedman era was stiffer than that of radiation, 
then the primordial backgrounds at high, low, and very 
low frequencies could be significantly stronger than O ~ 
3 x 10~ 14 , and could be strong enough to detect. 

9.2 Waves from Phase Transitions in the Early 
Universe 

A stochastic background could also have been produced 
by phase transitions in the early universe [104, 105|. No 
known phase transition would put its waves into the high- 
frequency band, and even hypothetical phase transitions, 
optimized at high-frequencies, can be only strong enough 
for marginal detection by advanced LIGO interferometers. 
The prospects for LISA are a little better: A strongly first- 
order electroweak phase transition could produce low - 
frequency waves strong enough for LISA to detect |l05|| . 

9. 3 Waves from Cosmic Strings 



If cosmic strings |106|[ were produced in the early universe 
in as large numbers as some theorists have suggested [ fL0[ 
|TT| , their vibrations would produce a gravitational wave 
spectrum that is frequency independent, fl g = const, from 
below the very low-frequency band where pulsar timing 
operates, through LISA's low-frequency band, and on into 
and through the high-frequency band. Theory suggests 



[107 that such waves could be as strong as fi s 
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level that is already being constrained by pulsar timing 
observations (Section 2.3). LIGO's first interferometers 
will operate at this same level, and by the time LIGO's 
advanced interferometers and LISA reach fl g ~ (a few) x 
10 -10 , pulsar timing might be in that same ballpark. 

To summarize: There are known mechanisms that could 
easily produce a measurable stochastic background in the 
high, low, and very-low frequency bands. However, the 
odds of the background being that large, based on cur- 
rently fashionable ideas, are not great. Despite this, a 
vigorous effort to detect background waves and to map 
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their spectrum will surely be made, since the cosmologi- 
cal implications of their discovery could be profound. 

§10 Conclusion 

It is now 35 years since Joseph Weber initiated his pioneer- 
ing development of gravitational-wave detectors f23] and 



25 years since Forward [10S] and Weiss [|T6| initiated work 
on interferometric detectors. Since then, hundreds of tal- 
ented experimental physicists have struggled to improve 
the sensitivities of these instruments. At last, success is 
in sight. If the source estimates described in this review 
article are approximately correct, then the planned inter- 
ferometers should detect the first waves in 2001 or several 
years thereafter, thereby opening up this rich new window 
onto the Universe. 
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